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INTRODUCTION 


The  computer  program  is  based  on  a  program  developed  by  Drs.  Walter 
Krauss  and  Peter  Sacher  at  Messerschmidt-Boelkow-Blohm  in  Munich, 
Germany,  and  reported  in  Reference  1.  The  MBB  Program  uses  a  method 
developed  by  Dr.  Paul  Rubbert  and  Gary  Saaris  at  the  Boeing  Company, 
Reference  2,  which  in  turn  stems  from  the  well-known  Douglas  Neumann 
Program  originated  by  John  Hess  and  A.  M.  0.  Smith,  Reference  3. 

A  listing  of  the  MBB  Program  was  provided  by  Dr.  Wolfgang  Schmidt  of 
the  Dornier  Company.  The  present  computer  program  retains  the  basic 
structure  of  the  MBB  Program,  but  has  been  extended  to  include  a 
plotting  package,  analysis  of  yawed  configurations,  and  many  other 
features  useful  in  the  analysis  of  bluff  bodies  in  subsonic  flow. 


1 


AERODYNAMIC  THEORY 


DESCRIPTION  OF  METHOD 


The  configuration  surface  is  divided  into  a  large  number  of  panels, 
each  of  which  contains  a  constant  source  distribution.  In  addition, 
an  internal  vortex  lattice  is  located  along  the  mean  chord  of  lifting 
surfaces  to  provide  circulation  to  the  flow.  A  typical  configuration 
panel  subdivision  is  shown  in  Figure  1. 

Analytical  expressions  for  the  perturbation  velocity  field  induced 
by  a  constant  source  distribution  on  an  arbitrary  quadrilateral  panel 
are  given  by  Hess  and  Smith  (Reference  3).  Similarly,  the  velocity 
field  induced  by  the  elements  of  a  vortex  lattice  are  given  by  Rubbert 
and  Saaris  (Reference  2).  The  perturbation  velocities  are  used  to 
calculate  the  coefficients  of  a  system  of  linear  equations  relating 
the  magnitude  of  the  normal  velocities  at  the  panel  control  points  to 
the  unknown  source  and  vortex  strengths.  The  source  and  vortex  strengths 
which  satisfy  the  boundary  condition  of  tangential  flow  at  the  control 
points  for  a  given  Mach  number  and  angle  of  attack  are  determined  by 
solving  this  system  of  equations  by  an  iterative  procedure.  The  pressure 
coefficients  at  panel  control  points  are  then  calculated  in  terms  of  the 
perturbation  velocity  components,  and  the  forces  and  moments  acting  on, 
the  configuration  obtained  by  numerical  integration. 

The  perturbation  velocity  components  induced  by  the  sources  and  vertices 
are  described  in  the  following  paragraphs,  together  with  the  formation 
and  solution  of  the  boundary  condition  equations,  and  the  procedure 
used  to  calculate  the  pressure  coefficients,  forces  and  moments  on  the 
configuration. 
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Figure  1.  Source  and  Vortex  Panel  Arrangement 
on  Wing- Body  Combination. 


THE  INCOMPRESSIBLE  VELOCITY  COMPONENTS 


The  perturbation  velocity  components  u,  v,  and  w  induced  by  a  constant 
distribution  of  sources  on  an  arbitrary  quadrilateral  panel  are  derived 
in  Reference  3,  so  only  the  final  expressions  will  be  reported  here. 
Consider  the  panel  shown  in  Figure  2. 


The  panel  is  assumed  to  lie  in  the  plane  z  ■  0,  and  the  corners  are 
numbered  clockwise  for  reference.  The  perturbation  velocities  cf  point 
P(x,y,z)  are  given  as  the  sum  of  the  contributions  of  the  four  sides  of 
the  quadrilateral  as  follows: 


u  * 

'H12Q12 

“H23Q23 

_H34Q34  "H41Q41 

(1) 

V  * 

K12Q12 

+K23Q23 

+K34^34  +K41«41 

(2) 

w  * 

f  1  «  -  J12 

_J23  "J34  'J41] 

(3) 

4 


where 


pij  "  (Ei  *  *>  Kij  +  (ni  -  yi  »y 

Tij  ■  (5i  - x)  Ku  +  ^  -  y>  Hy 

ri  "  "V  2+  -r>1>  2+z2i1/2 

*13  -  KCi-Ej)  2+  Oij-Mj)  2J1/2 

and  A0  *  2u  if  the  point  P  lies  inside  the  boundary  in  the  plane 
of  the  panel;  A0  -  0  otherwise. 

For  the  points  located  more  than  four  times  the  length  of  the  major 
diagonal  from  the  panel  centroid,  the  quadrilateral  is  approximated  by 
a  point  source  at  the  centroid.  This  simplifies  the  expression  for  the 
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velocity  components  considerably. 

In  this  case, 

u  *  (x  -  x)  S/r^  (M 

v  =  (y  -  y)  S/r^  (5) 

w  =  (z  -  z)  S/r^  (6) 

where 

r  =  [ (x  -  x)2  +  (y  -  y)2  +  (z  -  z)2]  1 ^2 

S  =  panel  area 

and  x,  y,  z  are  the  coordinates  of  the  panel  centroid. 

Additional  multipole  expansion  formulas  for  the  velocity  components 
given  in  Reference  3  are  not  used  in  this  program. 

The  perturbation  velocity  components  induced  by  a  line  vortex  are 
derived  in  Reference  2.  The  line  vortex  is  represented  by  a  vector 
as  shown  in  Figure  3  below.  ^It  induces  a  counterclockwise  circulation 
in  a  plane  perpendicular  to  L  if  the  vortex  strength  is  positive. 


Figure  3.  Line  Vortex  Geometry. 
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The  velocity  at  P(x,  y,  z,)  is  perpendicular  to  the  plane  containing 

# 

L  and  the  point  I  ,  and  is  given  by  Biot-Savart  s  Law  as 


4  7T 


f  £*t 

J  |dxl| 

Tx? 


sin  0 


D 


2 


ds 


4tt  of*  rL)2 


(cosQ.  -  cos  0 .  ) 


(7) 


A  quadrilateral  vortex  is  composed  of  four  vortex  segments  of  equal 
strength.  The  velocity  induced  by  a  quadrilateral  vortex  is 


u.i  +  v.j  +  w.k  ■  Vn 

n=] 


(8) 


where  V  is  the  velocity  induced  by  segment  n,  and  u,  v,  and  w  are 
n 

the  perturbation  velocity  components. 

A  vortex  lattice  consists  of  a  series  of  quadrilateral  vortices  of 
varying  strengths,  with  adjacent  edges  superimposed.  The  vortex  wake 
is  approximated  by  giving  the  last  vortex  quadrilateral  a  large  but 
finite  length.  The  net  strength  of  the  trailing  vortices  in  the  wake  is 
the  sum  of  the  strengths  of  the  individual  elements  in  the  lattice,  as 
indicated  in  Figure  4. 
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The  relative  strengths  of  the  individual  bound  vortices  in  the  lattice 
are  specified  in  advance.  The  circulation  around  each  airfoil  section 
is  determined  by  the  net  strength  of  each  vortex  lattice. 

COMPRESSIBILITY  CORRECTIONS 


The  velocity  components  in  compressible  flow  are  found  by  applying 
Gothert's  Rule  (Reference  A).  Two  options  are  available  in  the  program 
for  applying  the  compressibility  corrections,  and  are  designated  Rule  1 
and  Rule  2. 

Rule  1  applies  the  method  originally  proposed  by  Gothert.  The 
incompressible  velocity  components  are  calculated  on  an  analogous 
body  obtained  by  the  following  transformation: 

x  ■=  x 
a 

ya  =  By  (9) 


where 

B  =  ^1-M2 

The  boundary  conditions  of  tangential  flow  are  applied  on  the  analogous 
body,  and  the  resulting  incompressible  perturbation  velocities  are  trans¬ 


formed  back  to  the 

real  body  by 

u  * 

u  /  B2 
a 

v  * 

v  /  B 
a 

(10) 

w  = 

w  /  B 
a 

The  total  velocity 

vector 

at  a  given  point  is  then 

U  - 

U  cos 

OO 

a  cos  6  +  u 

V  = 

U  sin 
00 

B  +  v 

(ID 

w  = 

U  sin 
00 

a  cos  B  +  w 

it  is  now  known  that  this  compressibility  rule  yields  good  results 
only  for  slender  bodies  at  small  angles  of  attack.  The  validity  of 
this  rule  decreases  with  increasing  values  of  the  surface  slope.  This 
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effect  is  particularly  noticeable  for  two-dimensional  airfoil  sections. 

In  the  vicinity  of  the  nose,  Gothert's  Rule  (which  is  equivalent  to  the 
Prandtl-Glauert  Rule  in  this  example)  gives  excessively  high  suction 
peaks  on  the  upper  surface.  The  reason  for  this  failure  of  the  theory 
is  the  manner  in  which  the  boundary  conditions  are  satisfied.  Since 
the  boundary  conditions  are  satisfied  at  the  surface  of  the  analogous 
body  which  is  thinner  by  the  factor  B  than  the  real  body,  the  curvature 
of  the  flow  near  the  nose  is  correspondingly  increased,  resulting  in 
higher  suction  peaks.  In  order  to  eliminate  this  effect,  it  is  necessary 
to  satisfy  the  boundary  conditions  on  the  surface  of  the  real  body. 

Rule  2  was  first  proposed  by  Kraus  in  Reference  1.  Beginning  with 
the  analog  body  as  before,  the  expressions  for  the  perturbation  velocity 
components  are  corrected  for  compressibility,  using  Equation  10,  prior 
to  solving  the  boundary  condition  equations.  The  boundary  conditions  of 
tangential  flow  are  then  applied  on  the  surface  of  the  real  body, 
resulting  in  improved  results  for  the  velocities  and  pressure 
coefficients. 

THE  BOUNDARY  CONDITION  EQUATIONS 

The  boundary  condition  of  tangential  flow  at  panel  control  points 
establishes  a  system  of  linear  equations  for  determining  the  strengths 
of  the  source  and  vortex  distributions.  The  geometrical  relationship 
between  each  panel  and  control  point  is  required  to  evaluate  the 
coefficients  of  this  system  of  equations. 

Panel  Geometry 

A  typical  panel  subdivision  of  a  wing-body  configuration  is  illustrated 
in  Figure  1.  A  reference  coordinate  system  is  established  with  origin 
at  or  near  the  nose  of  the  configuration,  having  an  x-axis  lying  in  the 
plane  of  symmetry  parallel  to  the  body  axis,  and  a  vertical  z  axis. 
Symmetry  of  the  body  about  the  x,  z  plane  is  not  required.  However, 
if  the  body  is  symmetric,  only  those  panels  located  on  the  positive 
y  side  of  the  x,  z  plane  are  required. 

The  body  panel  corners  are  defined  by  the  intersections  of  a  series  of 
planes  normal  to  the  x  axis,  and  longitudinal  meridian  lines.  A 
maximum  of  70  body  sections  may  be  used,  and  each  section  may  contain 
up  to  60  points  around  the  half-circumference.  The  body  panel  corner 
points  may  be  shifted  longitudinally  to  aid  in  panelirg  wing-body 
intersections.  The  body  panels  are  numbered  in  sequence  from  the  top 
to  the  bottom  of  each  circumferential  ring,  starting  from  the  most 
forward  ring. 
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The  wing  panels  are  defined  by  the  intersections  of  a  series  of  vertical 
planes  parallel  to  the  plane  of  symmetry,  and  lines  of  constant  percent 
chord.  A  maximum  of  40  wing  sections  may  be  defined,  each  containing 
up  to  60  airfoil  ordinates.  The  same  number  of  ordinates  are  required 
on  the  upper  and  lower  surfaces  of  the  airfoil,  at  approximately  the 
same  percent  chord  locations,  in  order  to  properly  define  the  internal 
vortex  panels. 

The  wing  panel  corner  points  may  be  shifted  laterally  to  aid  in 
paneling  wing-body  intersections.  The  wing  panels  are  numbered 
sequentially,  and  follow  the  body  panels.  Beginning  with  the  inner 
chordwise  strip  of  panels,  the  numbering  starts  at  the  trailing  edge 
of  the  lower  surface,  and  ends  at  the  trailing  edge  of  the  upper 
surface.  A  maximum  of  1500  panels  may  be  used  to  define  the  external 
surfaces  of  the  wing  and  body.  Ii  the  configuration  is  symmetric, 
this  implies  a  maximum  of  750  panels  on  one  side  of  the  x,  7.  plane. 

Vortex  lattice  panels  are  automatically  defined  on  the  mean  chord 
plane  of  the  wing.  The  panel  corner  points  are  obtained  by  averaging 
the  upper  and  lower  surface  airfoil  ordinates  at  each  percent  chord 
location.  One  additional  vortex  panel  is  defined  in  the  wake  aft  of 
the  trailing  edge  of  each  chordwise  strip  of  wing  panels  to  provide 
control  points  for  satisfying  the  Kutta  condition.  The  additional 
panel  lies  in  the  plane  of  the  trailing  edge  bisector.  For  wing-body 
combinations,  additional  vortex  lattices  are  required  inside  the  body 
to  provide  a  mechanism  for  carry-over  of  lift.  A  maximum  of  35  vortex 
lattices  may  be  defined,  and  these  are  numbered  subsequently  following 
the  wing  panels. 

The  four  input  points  defining  a  panel  do  not  necessarily  lie  in  the 
same  plane.  The  technique  used  to  approximate  the  panel  by  an  equiva¬ 
lent  planar  panel  was  developed  by  Hess  and  Smith,  Reference  5,  and  is 
summarized  in  Appendix  I.  Using  this  method,  a  panel  coordinate  system 
is  defined  with  origin  at  the  panel  centroid  and  lying  in  the  mean 
plane  of  the  input  points.  The  x  axis  of  the  panel  coordinate  system 
is  parallel  to  one  diagonal,  the  z  axis  is  normal  to  the  plane  of  the 
panel,  and  the  y  axis  is  perpendicular  to  the  other  two.  Since  the 
velocity  components  induced  by  the  source  distributions  are  given  in 
terms  of  the  panel  coordinate  system,  a  nine  element  transformation 
matrix  Tij  is  calculated  for  each  panel  to  transform  the  coordinates 
of  points  and  the  components  of  vectors  from  the  reference  coordinate 
system  to  the  panel  coordinate  system.  In  addition,  the  panel  area, 
the  coordinates  of  the  centroid,  and  the  length  of  the  principal  diagonal 
are  calculated. 
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Normal  Velocity  at  Panel  Control  Points 


Each  surface  panel  is  assigned  a  control  point  located  at  the  panel 
centroid.  Each  vortex  lattice  is  assigned  a  control  point  just  behind 
the  trailing  edge  of  the  wing  in  the  plane  of  the  trailing  edge 
bisector.  (This  point  is  normally  located  1  percent  of  the  local  chord 
behind  the  trailing  edge.) 

The  resultant  velocity  normal  to  panel  i  at  its  control  point  is  the 
sum  of  the  normal  component  of  the  free-stream  velocity  vector  and  the 
normal  velocities  induced  by  the  panel  source  and  vortex  distributions. 
Setting  the  magnitude  of  the  free-stream  velocity  vector  equal  to  unity, 
its  component  normal  to  panel  i  is 

R.  =  cos  a  •  cos  8  •  n  +  sin  8  •  ny,  +  sin  a  •  cos  8  •  n  (12) 

xi  Zi 

where  n  ,  n  ,  and  n  are  the  direction  cosines  of  the  normal  of 
Xi  yi  Zi 

panel  i  (see  Appendix  I),  a  is  the  angle  of  attack  and  8  is  the  angle 
of  yaw  of  the  free-stream  velocity  vector  in  the  reference  axis  system. 


The  normal  component  of  velocity  induced  at  the  control  point  of  panel 
i  by  the  source  and  vortex  distributions  is  given  by 

N 


->: 

j=i 


(n 


+  n 


)  o, 


ij 


(13) 


'ij 


where  v  ,  v  ,  and  v  are  the  three  components  of  velocity  parallel 
Xij  yij  Zij 

to  the  reference  axis  at  control  point  i  induced  by  a  unit  strength 
source  or  vortex  distribution  on  panel  j  and  Oj  is  the  strength  of  the 

singularity. 


The  three  components  of  velocity  parallel  to  the  reference  axes  are 
obtained  by  multiplying  the  velocity  components  given  by  Equations 
(4),  (5),  and  (6)  in  the  panel  coordinate  system  by  the  transformation 
matrix  given  in  Appendix  I.  For  example, 


ij 

Uij 

Xi 

+ 

vij 

tlyij 

+ 

Wij 

‘lzij 

ij 

*  Uij 

t2*il 

+ 

vij 

t2yij 

+ 

Wij 

(14) 
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xij 

ij 

yij 

ij 

Zij 
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Combining  Equations  (12)  and  (13), 


n, 


Ri  +  Ai 


r±+)2  a-  a 


J-i 


ij  j 


(15) 


where  the  aerodynamic  influence  coefficient  a^  is  given  by  Equation (13) 
Solution  of  the  Boundary  Condition  Equations 

The  boundary  condition  of  tangential  flow  at  panel  control  points  is 
satisfied  if  the  normal  velocities  are  set  equal  to  zero  on  all 
panels. 


Thus 


or 


N 

1  Vni  '  ° 

i=l 

N  N 

2  I  •«  '  ■  2  Ri 


L* 

i-1  J-l 


N 

y 

j*i 


(16) 


In  matrix  notation. 


tA±J  ]  {Oj}  =-{Ri) 


(17) 


where  is  the  matrix  of  aerodynamic  influence  coefficients,  and 
the  right  side  of  the  equation  is  given  by  Equation  (12). 


This  system  of  equations  can  be  solved  by  direct  inversion  to  determine 
the  unknown  source  and  vortex  strengths.  However,  for  the  large  order 
matrices  usually  encountered  in  aerodynamic  problems,  an  iterative 
solution  procedure  is  given  in  Reference  3.  A  modified  Gauss-Siedel 
iteration  procedure  is  employed  in  this  computer  program. 
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This  matrix  is  subdivided  into  four  partitions  as  follows: 


where  A  is  the  matrix  gi'r»->3  the  influence  of  the  source  panels  on 
ss 

the  surface  control  points. 

Agv  is  the  matrix  giving  the  influence  of  the  source  panels  on  the 
vortex  lattice  control  points. 

A  is  the  matrix  giving  the  influence  of  the  vortex  lattices  on  the 
vs 

surface  control  points,  and 

A^  is  the  matrix  giving  the  influence  of  the  vortex  lattices  on  the 
vortex  lattice  control  points. 

Equation  (17)  may  now  be  written  as 


A 

ss 

Avs 

• 

Si- 

m 

R  . 
si 

(18) 

A 

sv 

A 

w 

*i 

R  . 
vi 

A  s . 
ss  i 

+  A  y. 
vs  'i 

"  -Rsi 

(19) 

A  s . 
sv  i 

+  A  y. 
w  i 

=  -  R 

Vi 

(20) 
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where  are  the  unknown  source  strengths,  and  y^  are  the  unknown 

vortex  strengths,  Agg  is  a  square  matrix  of  order  NS,  and  is 

a  square  matrix  of  order  NL,  where  NL  is  generally  much  smaller  than 
NS. 


The  first  step  in  each  of  the  iteration  cycles  is  to  use  Equation 
(19)  only.  The  values  for  y^  are  taken  from  the  previous  cycle  (or 


set  equal  to  zero  on  the  first  cycle)  and  a  solution  for  the  array 
{s^}  is  obtained  by  the  Gauss-Seidel  procedure.  These  values  for  s^ 
are  then  used  in  Equation  (20)  to  obtain  y^  by  direct  inversion: 


y 


i 


(21) 


These  values  are  now  used  in  the  first  step  of  the  next  cycle,  and 
the  procedure  continues  until  convergence  is  achieved.  The  criterion 
for  convergence  is 


SS  •,  ^  * 

i-l\J-l  1J  £ 1 


VSij 


+  R 


<  e 


(22) 


where  e  is  some  small  number  specified  by  the  user.  Normally  e  <  10 


-3 


More  elaborate  iteration  schemes  using  smaller  partitions  of  the  A 

s  s 

matrix  are  <  ?scribed  in  Reference  6,  but  have  not  been  incorporated 
into  the  present  program. 


CALCULATION  OF  THE  PRESSURES,  FORCES,  AND  MOMENTS 

Once  the  source  and  vortex  strengths  have  been  determined,  the  three 
components  of  velocity  at  control  point  i  may  be  obtained. 


N 


cos  a  cos  3+/  v  a. 

j4  -13  1 

N 

(23) 

sin  8  +  ;  v  0. 

^  ?ii  J 

J-l  J 

(24) 
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sin  ot  cos 


(25) 


w. 


N 

b  +n: 


& 


V  0. 

zij  j 


where  the  0,  includes  both  source  and  vortex  strengths,  and  v  , 
j  Xij 

v  ,  and  v  are  defined  following  Equation  (13).  The  pressure 

Yij  Zij 


coefficient  is  calculated  using  the  exact  isentropic  formula 


2 

=  —  -  - 
w2 

1  K-l  w2  ,,  2. 

1  +  — —  M  (l-q±) 

3.5 

-  1 
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(26) 


where 
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qi 


2 

U1 


+ 


w 


2 

i 


For  M  <  .1, 


the  program  uses  the  simpler  formula 


(27) 


The  forces  and  moments  acting  on  the  configuration  can  now  be  obtained 
by  numerical  integration.  The  normal  force,  side  force,  axial  force, 
and  pitching  moments  (about  the  origin  of  coordinates)  of  panel  i  are 
given  by 


Xi  " 

Si  CPi 

n 

Xi 

(28) 

Yi  * 

Si  CPl 
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yi 

(29) 

Zi  - 

Si  CPi 
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ZiYi  - 

Vi 
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where  S.  is  the  area  of  the  panel,  n  ,  n  ,  and  n  are  the  direction 

1  xi  y± 

cosines  of  the  normal,  and  x^,  and  are  the  coordinates  of  the 

panel  control  point. 

The  total  force  and  moment  coefficients  are  obtained  by  summing  the 
panel  forces  and  moments  on  both  sides  of  the  plane  of  symmetry 


C 


Z 


'"Jr* 

H. 
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i-1 
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1 

S  c 
w 
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i-1 
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S  c 
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(34) 


(35) 


(36) 


(37) 


(38) 


(39) 


Finally,  the  lift,  side  force, and  drag  coefficients  are 


CL 

*  C^cos  a  - 

(C^cosB  -  C^sinS)  sin  a 

(40) 

cs 

■  C^cos  8  + 

C^sin  8 

(41) 

CD 

-  (C^cos  8  - 

■  CYsin8)  cos  a  +  C^sinoi 

(42) 
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The  program  computes  the  forces  and  moments  acting  on  the  body  and 
the  wing,  and  sums  them  to  obtain  the  total  forces  and  moments  of  the 
configuration.  In  addition,  wing  section  forces  and  moments  may  be 
calculated  at  the  user's  option. 

SEPARATED  FLOW  MODEL 


The  flow  external  to  the  boundary  layer  and  the  separated  wake  is 
essentially  potential  flow.  To  obtain  a  mathematical  model  with 
potential  flow  everywhere,  the  boundary  layer  and  separated  wake 
in  the  real  case  are  replaced  by  fluid  originating  from  the  body 
surface  as  shown  in  Figure  5.  Rules  for  the  distribution  of  surface 
normal  velocity  to  account  for  boundary  layer  growth  upstream  of 
separation  have  been  formulated,  and  their  use  requires  matching  the 
boundary  layer  and  the  potential  flow  solution  by  iteration.  Rules 
for  distributing  surface  normal  velocity  in  the  separated  region  to 
obtain  fluid  which  will  displace  the  potential  flow  originating  from 
upstream  in  the  same  way  as  the  separated  wake  in  the  real  case  have 
not  been  formulated.  However,  for  very  unstreamlined  bodies,  a  plaus¬ 
ible  approach  is  to  make  the  surface  normal  velocity  equal  to  the  free- 
stream  velocity  component  normal  to  the  surface: 


vs  *  n  =  •  n  (43) 

Vg  *  surface  velocity 

-V 

=  free-stream  velocity 

— ► 

n  =  unit  surface  normal 

To  calculate  the  location  of  separation  requires  a  matching  of  the 
boundary  layer  and  potential  flow  solutions  similar  to  that  used  for 
the  boundary  layer  growth,  a  formidible  job.  However,  for  an  aft 
body  shape  with  rapid  closure  or  for  bodies  having  sharp  corners, 
making  an  intuitive  estimate  of  the  separation  point  is  consistent 
with  the  use  of  Equation  43. 
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SEPARATED  — 
WAKE  J 


& 


BOUNDARY  LAYER 


(a)  REAL  FLOW 


Figure  5.  Modeling  of  Potential  Flow  to  Account  for  Boundary 

Layer  and  Wake. 
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This  approximate  separation  modeling  was  applied  to  the  BO  105  helicopter 
fuselage  configuration.  Figure  6  shows  the  different  choices  used  for 
the  region  where  fluid  was  ejected  according  to  Equation  A3.  The  choice 
labeled  Case  A  gave  the  best  agreement  with  experiment  upstream  of 
separation  as  shown  in  Figures  7  and  8.  The  pressure  distributions  for 
Waterline  10  (see  Figure  8)  show  that  the  ejection  region  was  extended 
too  far  forward  for  Case  B.  The  top  centerline  data  (see  Figure  7) 
shows  that  for  Case  B  the  ejection  region  was  extended  too  high  on 
the  body.  The  presence  of  the  boom  and  a  milder  closure  evidently 
eliminates  separation  on  the  upper  part  of  the  body. 

It  should  be  noted  that  by  proper  choice  of  outflow  distribution,  the 
potential  flow  model  shown  in  Figure  5b  can  make  the  flow  external  to 
the  dividing  streamline  agree  well  with  the  real  case  and  hence  give 
the  proper  surface  pressure  upstream  of  separation,  but  it  cannot  be 
expected  to  simultaneously  provide  surface  pressures  that  agree  with 
the  real  flow  in  cr.e  separated  region.  The  irrelevant  result  in  this 
region  is  illustrated  in  Figure  8  beyond  X  ■  39.  An  approximate  replace¬ 
ment  is  to  use  the  pressure  at  separation  over  the  entire  separated  sur¬ 
face  region.  The  experimental  results  shown  in  Figure  8  suggest  this 
approach.  This  level  can  be  estimated  from  the  potential  flow 
solution  using  boundary  layer  separation  criteria.  Applying  a  modified 
Townsend  criterion  to  the  Waterline  10  pressure  distribution  for 
separation  modeling  Case  A  gives  the  pressure  level  shown  as  C  in 

*s 

Figure  8.  This  modified  Townsend  criterion,  developed  by  F.  A.  Dvorak, 
is; 


C 

P 

s 


-  pressure  coefficient  at  separation 

Re 
x 

.  o 


U-C  ) 


-83.961  +  38.645  log 


(1-C  ) 


+  C 
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Re 


=  pressure  coefficient  at  recovery 

=  boundary  layer  development  length  to  recovery  point 
*  recovery  pressure  gradient 
=  Reynolds  number  based  on  length  xq 


skin  friction  coefficient  at  recovery 
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x  values  of  panel  edges 


Waterline  10  (z«4) 


Regions  of  ejection  for  modeling  separation: 


Region  for  Case  A  is  shown  by  □ 

Region  for  Case  B  includes  that 
for  Case  A  and  region  shown  by 

Boundary  Condition  enforced  on  ejection  regions: 


Figure  6.  The  BO  105  Helicopter  Fuselage  Showing  Paneling 
and  Separation  Modeling. 
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Figure  7.  Pressure  Distribution  Along  Top  Centerline  of  the  BO  105. 


SEPARATION  PRESSURE  COEFFICIENT,  Cps,  WAS  CALCULATED  FROM  CASE  A  PRESSURE 
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Figure  8.  Pressure  Distribution  Along  Waterline  10  of  the  BO  105. 


The  values  for  C  ,  C  ,  and  x  that  were  taken  from  the  Case  A  pressure 
p  p  o 

*o  r 

distribution  to  calculate  C  are  indicated  in  Figure  8.  The  boundary 

Ps 

layer  development  length  x  was  approximated  by  using  the  axial  dis¬ 
tance  from  the  nose  of  the  body  to  the  recovery  point.  For  the  recovery 
pressure  gradient  C',  the  maximum  gradient  in  the  recovery  region  was 
used.  P 

For  the  cases  shown  in  Figures  7  and  8,  experimental  results  are  avail¬ 
able  to  determine  the  best  ejection  configuration  for  the  paneling  used, 
namely,  Case  A.  An  obvious  problem  is  that  in  the  absence  of  experi¬ 
mental  data,  it  is  not  known  where  to  start  the  surface  ejection.  It 
may  be  possible  to  determine  this  point  by  an  iterative  procedure  using 
a  boundary  layer  separation  criterion.  The  idea  is  to  start  the 
surface  ejection  at  the  separation  point  calculated  from  the  preceding 
cycle,  an  initial  guess  being  used  for  the  first  cycle.  The  example 
in  Figures  7  and  8  shows  that  for  a  jump  start,  ejection  must  begin 
some  distance  downstream  from  the  calculated  separation  point.  Clearly, 
the  iterative  procedure  suggested  requires  a  gradual  initiation  of 
ejection,  and  the  shape  used  for  this  initiating  "ramp"  is  critical  for 
convergence  and  for  obtaining  an  accurate  modeling  of  separation  as 
pictured  in  Figure  5.  Also  denser  paneling  than  shown  in  Figure  6  is 
necessary  to  provide  more  resolution  for  the  ejection  distribution. 
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COMPUTER  PROGRAM 


PROGRAM  DESCRIPTION 


The  program  developed  to  calculate  the  pressure  distribution  and 
aerodynamic  characteristics  of  wing-body  combinations  in  subsonic  flow 
is  written  in  FORTRAN  IV.  A  maximum  of  1500  source  panels  and  35  vortex 
lattices  may  be  used  to  represent  the  configuration.  It  is  designed  to 
operate  on  both  the  CDC  6600  or  IBM  360/370  series  of  computer  with 
minor  modifications.  The  program  requires  approximately  210,000 
(octal)  words  storage  on  the  CDC  computer,  and  operates  in  OVERLAY 
mode.  The  program  requires  five  peripheral  disc  files  in  addition 
to  the  input  and  output  files. 

PROGRAM  STRUCTURE 


The  overlay  structure  of  the  program  is  illustrated  in  Figure  9.  The 
main  overlay  program  is  designated  WBOLAY,  and  calls  the  three  primary 
overlay  programs  WBPAN,  WBPLOT,  dnd  UBAEOR.  The  complete  program  consists 
of  16  subroutines  in  addition  to  standard  library  and  plot  subroutines. 
Descriptions  of  these  subroutines  are  contained  in  Appendix  II  of  this 
report . 

PROGRAM  INPUT  DATA 


The  input  to  the  program  is  divided  into  three  parts:  the  geometry 
input,  the  plot  input,  and  the  aerodynamic  input.  The  input  require¬ 
ments  of  each  part  are  described  below.  A  sample  input  is  given  in 
Appendix  III,  and  a  sample  output  in  Appendix  IV. 


Figure  9.  Program  Overlay  Structure. 
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Geometry  Input  Cards 


If  the  configuration  is  symmetrical  about  the  x,  z  plane,  geometrical 
input  is  required  for  only  one  side  of  the  configuration.  The  convention 
used  herein  is  to  present  that  half  of  the  configuration  lying  on  the 
positive  y  side  of  the  x,  z  plane.  If  the  configuration  is  not  symmetric, 
complete  geometrical  input  is  required. 

Card  1  -  General  Identification  -  Card  1  contains  any  desired  identifying 
information  in  Columns  1-80. 


Card  2  -  Configuration  Parameters  - 


Col.  Variable  Value 

10  CASE  1 

2 

3 


Description 

Isolated  body  only 
Isolated  wing  only 
Wing-body  combinations 


20  PLOT  0 

1 


No  plot  output 
Plot  output  requested 


30  SIM  0 

1 

-1 


Configuration  symmetric  about  x,  z 
plane;  panel  geometry  required  on  one 
side  only 

Configuration  symmetric  about  x,  z 
plane.  Panel  geometry  input  required 
on  one  side  only;  panel  geometry  output 
calculated  for  both  sides.  (Used  when 
analyzing  symmetric  configuration  in 
yaw.) 

Unsymmetric  configuration.  Panel 
geometry  input  required  for  both  sides 


40  ISAVE  0 

1 

-1 


Geometry  and  influence  coefficient 
matrices  not  saved 
Geometry  and  influence  coefficient 
matrices  saved  in  previous  run  to  be 
used  (TAPE  11  must  be  requested  and 
card  sets  3-8  omitted) 

Geometry  and  influence  coefficient 
matrices  to  be  saved  on  TAPE  11 
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Col.  Variable  Value 


Description 


50  PRINT  0 

1 

2 

3 

4 


Normal  output  -  see  Program  Output 
Section  (page  41  ) 

Optional  output  1  -  Includes  panel 
geometry,  coordinate  transformation 
matrices,  and  panel  forces  and  moments 
Optional  output  2  -  Panel  velocity 
components  and  influence  coefficients 
Requires  large  line  count  limit 
Optional  output  3  -  The  aerodynamic 
influence  coefficient  matrix,  the  right 
side  of  the  matrix  equation,  and  all 
solution  iterations 

Optional  output  4  -  This  option  prints 
out  the  successive  solution  iterations 
only 


Note:  The  normal  output  is  always  printed  in  addition  to  any  optional 

output  selected. 

Card  Set  3  -  Single  Panel  Input  -  This  card  set  allows  individual  panels 
to  be  input  by  specifying  the  coordinates  of  the  four  corner  points  in 
clockwise  order.  Any  number  of  panels  may  be  input  in  this  manner.  It 
also  allows  individual  panels  to  be  deleted  by  specifying  the  panel 
indices.  A  maximum  of  100  panels  may  be  deleted. 


Card  3A 

-  Single 

Panel  Control  Card 

Col. 

Variable 

Value 

Description 

1-10 

SINGPA 

0 

1 

No  single  panel  input;  omit  card  set 
3B,  continue  reading  input  cards 

Corner  point  coordinates  of  this  panel 
follow  on  card  set  3B 

11-20 

NOPAN 

Arbitrary 

Integer 

Number  of  panels  to  be  deleted.  If 
non-zero,  panel  indices  follow  on  card 
set  3C 

Card  3B 

-  Panel  Corner  Point  Input 

Col. 

Variable 

Value 

Description 

1-10 

X(I) 

Arbitrary 

(floating 

point) 

x  coordinate  of  corner  I 

11-20 

Y(I) 

II 

y  coordinate  of  corner  I 
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Col. 

Variable 

Value 

Description 

21-30 

Z(I) 

Arbitrary 

z  coordinate  of  corner  I 

(floating 

point) 

Repeat  card  3B  four  times,  once  for  each  corner  of  the  panel. 
Card  Set  3C  -  Indices  of  Deleted  Panels 


NOPAN  indices  of  deleted  panels  are  read  (7110)  format)  if  NOPAN  >  0 
on  card  3A.  A  maximum  of  100  panels  may  be  deleted.  Wing  and  body 
vortex  lattice  control  panels  may  not  be  deleted. 

Card  Set  4  -  Body  Panel  Input  -  This  card  set  allows  the  body  panels 
to  be  calculated  automatically,  from  the  section  geometry  data.  Five 
options  are  available  for  inputting  the  section  geometry.  The  XYZ 
program  input  referred  to  below  conforms  with  the  format  of  Reference 
7.  Omit  this  card  set  if  CASE  *  2  on  card  2. 


Card  4A 

-  Number  of  Body  Sections 

Col. 

Variable 

Value 

Description 

1-10 

NB 

Arbitrary 

Integer 

Number  of  body  sections 
(2  <  NB  ;<  70) 

Card  4B 

-  Body  Section  Geometry 

Col. 

Variable 

Value 

Description 

1-10 

XBE 

Arbitrary 

(floating 

point) 

x  coordinate  of  origin  of  body  section  co¬ 
ordinate  system  except  blank  when  XYZ  pro¬ 
gram  input  format  is  used 

11-20 

YBE 

II 

Similarly  the  y  coordinate 

21-30 

ZBE 

II 

Similarly  the  z  coordinate 

31-40 

MB 

Arbitrary 

Integer 

No.  of  input  points  on  section 
(3  <  MB  <  60) 

If  MB  <0,  XYZ  program  input 
format  requested 

50 

OPT 

0 

Body  section  geometry  input  by  y-z 
coordinates  on  card  set  4C  and  4D 

1 

Body  section  geometry  same  as 

preceding  body  section  -  card  set 
4C  or  4D  omitted.  Note:  YBE  and  ZBE 
are  additive  to  preceding  values* 
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Card  4B  -  Body  Section  Geometry  (cont’d) 


Col. 

Variable 

Value 

Description 

2 

Body  section  geometry  input  in 
polar  coordinates  r,  0  on  card  set  4C 

3 

Body  of  revolution,  section  geometry 
input  as  section  radius  and  theta 
increment  on  card  set  4C 

60 

FLAG 

0 

Normal  body  section 

1 

Terminal  body  section  (end  of  current 
body  panel  network) 

Card 

Set  4C  -  Body 

Section  Coordinates  (Normal  Input) 

Col. 

Variable 

Value 

Description 

1-10 

B(J) 

Arbitrary 

(floating 

point) 

y  coordinate  of  point  J  if  OPT  ■  0, 
or,  angular  coordinate  (in  degrees)  of 

J  if  OPT  -  2  or,  increment  angle  A9 
in  degrees  if  OPT  *  3 

11-20 

A(J) 

If 

z  coordinate  of  point  J  if  OPT  -  0, 
or,  r  coordinate  of  point  J  if  OPT  ■  2, 
or,  body  section  radius  if  OPT  =  3 

21-30 

D(J) 

II 

A  x  shift  of  point  J  if  OPT  ■  0  or  2 

Card  set  4C  contains  MB  cards  if  OPT  =  0  or  2,  contains  only  1  card  if 
OPT  -  3,  and  is  omitted  if  OPT  ■  1  or  MB  <0  on  card  4B.  See  Figure  11. 

Card  Set  4D  -  Alternate  XYZ  Input 


Col. 

Variable 

Value 

Description 

1-12 

D(J) 

Arbitrary 

x  coordinate 

of 

section 

13-24 

B(  J) 

(floating  point) 

It 

y  coordinate 

of 

point 

J 

25-36 

A(J) 

II 

z  coordinate 

of 

point 

J 

This  card  set  is  omitted  unless  MB  <0  In  Card  4B. 


Note:  Repeat  card  4B  and  card  sets  4C  or  4D  NB  times  to  complete 
card  set  4. 
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Card  Set  5  -  Wing  Panel  Input  -  This  card  set  allows  the  wing  and  vortex 
lattice  panels  to  be  calculated  automatically  from  the  wing  section 
data.  Three  options  are  available  for  inputting  the  wing  section 
geometry.  The  XYZ  program  input  referred  to  below  conforms  with  the 
format  of  Reference  7.  Omit  this  card  set  if  CASE  *  1  on  card  2. 


Card  5A 

-  Number 

of  Wing  Sections 

Col. 

Variable 

Value 

Description 

1-10 

NV 

Arbitrary 

Integer 

Number  of  wing  sections  (2<  NW  <  40) 

20 

KOORD 

1 

2 

Wing  section  ordinates  input  in 
percent  of  local  chord 

Wing  section  ordinates  input  are 
not  normalized 

Card  5R 

-  Wing  Section  Geometry 

Col. 

Variable 

Value 

Description 

1-10 

XBE 

Arbitrary 
(floating  point) 

x  coordinate  of  origin  of  wing  section  co¬ 
ordinate  system  except  blank  when  XYZ  pro¬ 
gram  input  format  is  used 

11-20 

YBE 

If 

Similarly  the  y  coordinate 

21-30 

ZBE 

ft 

Similarly  the  z  coordinate 

31-40 

CHRD 

If 

Chord  length  of  section 

41-50 

ALF 

II 

Section  twist  angle  (degrees)  (twist 
positive  for  dz/dx  negative  ) 

51-60 

XAL 

II 

Center  of  twist  in  percent  chord 

61-65 

MW 

Arbitrary 

Integer 

Number  of  coordinates  in  section  (5  <_  MW 
^  59)  .  Always  odd  number  if  internal  vor¬ 
tices  selected  If  MW  <  0,  XYZ  input  format 
requested. 

70 

OPT 

0 

Wing  section  ordinates  to  be  used  from 

card  set  5D  and  5E 


1  Wing  section  ordinates  same  as  preced¬ 

ing  section  -  card  set  5D  and  5E 
omitted 

75  FLAG  0  Normal  case  -  vortex  lattice  panels 

calculated  automatically 

Terminal  wing  section  (end  of  current 
wing  panel  network) 

30 
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Card  5B-Wing  Section  Geometry  (cont'd) 


Col.  Variable  Value 

2 

3 

80  DEL  0 

1 

Note:  The  values  of  CHRD,  ALF,  and 
KOORD  =  1  on  Card  5A. 

Card  Set  5C  -  Wing  Dihedral  Input 

Col.  Variable  Value 

I- 10  DELTA  Arbitrary 

(floating  point) 

II- 20  YO  " 

21-30  ZO 


Description 

No  vortex  lattice  panels  calculated 
for  this  section 

The  coordinates  of  the  last  bound 
vortex  in  the  vortex  lattice  are  read 
in  on  card  7  for  this  section 

No  wing  dihedral 

Dihedral  input  on  card  set  5C 

XA1  on  card  5B  are  required  only  if 

Description 

Dihedral  angle  (degrees) 

y  and  z  coordinates  of  axis  of 
rotation  of  wing  panel 


Omit  card  set  5C  if  DEL  «  0  in  card  5B 

Card  Set  5D  -  Wing  Section  Coordinates  (Normal  Input) 


Col. 

Variable 

Value 

Description 

1-10 

B(J) 

Arbitrary 

(floating 

x  coordinate  of  point  J 

point) 

11-20 

A(J) 

tl 

z  coordinate  of  point  J 

21-30 

C(J) 

11 

Vortex  lattice  strength  at  point  J 

31-40 

D(J) 

If 

A  y  shift  of  point  J 

Card  set  5D  contains  MW  cards  if  OPT  =  0,  and  is  omitted  if  OPT  «*  1 
or  MW<0  on  card  5B.  See  Figure  12. 
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Card  Set  5E  -  Alternate  XYZ  Input 


Col. 

Variable 

Value 

Description 

1-12 

B(J) 

Arbitrary 
(floating  point) 

x  coordinate  of  point  J 

13-24 

D(J) 

II 

y  coordinate  of  point  J 

25-36 

A(J) 

II 

z  coordinate  of  point  J 

51-60 

C(J) 

If 

Vortex  lattice  strength  at  point  J 

This  card  set  is  omitted  unless  MW<0  on  card  5B. 

Note:  Repeat  card  5B  and  card  sets  5C,  5D,  or  5E  NW  times  to  complete 
card  set  5. 

Card  Set  6  -  Vortex  Lattice  Control  Point  Location 


Col. 

Variable 

Value 

Description 

1-10 

WAKE 

Arbitrary 

Extension  of  vortex  lattice  into 

(floating  point)  wake  in  percent  chord 

11-20 

POINT 

II 

Location  of  vortex  lattice  control 

point  in  percent  chord  behind  trail¬ 
ing  edge 


Note:  These  values  are  not  used  if  FLAG  =  3  on  Card  5B.  Omit  card  6 
if  CASE  »  1  on  Card  2. 

Card  Set  7  -  Relocation  of  Vortex  Lattice  Terminal  Points  -  This  card 
set  is  omitted  unless  FLAG  =  3  on  card  5B.  For  each  wing  section  having 
FLAG  =  3,  two  additional  cards  are  required  to  specify  the  terminal  points 
of  the  streamwise  vortices. 

Card  Set  7A  -  Inboard  Terminal  Points 


Col. 

Variable 

Value 

Description 

1-10 

XLP 

Arbitrary  x  coordinate  of  inboard 

(floating  point)  lattice  terminal  point 

edge  of 

11-20 

YLP 

II 

y  coordinate  of  inboard 
lattice  terminal  point 

edge  of 

21-30 

ZLP 

II 

z  coordinate  of  inboard 
lattice  terminal  point 

edge  of 

32 


Card  Set  7B  -  Outboard  Terminal  Points  -  Same  as  card  7A  for  outboard 
edge  of  lattice  terminal  point. 

Card  Set  8  -  Body  Vortex  Lattice  Input  -  This  card  set  allows  additional 
vortex  lattices  to  be  located  inside  the  body  of  wing-body  combinations, 
and  is  omitted  if  CASE  <3  on  card  2. 

Card  Set  8A-  Number  of  Streamwise  Vortices  in  Body  Vortex  Lattice 
Network 


Col. 

Variable 

Value 

Description 

1-10 

NV 

Arbitrary 

Integer 

Number  of  streamwise  vortices 
in  body  vortex  lattice  network 
(NV  <  40) 

Note: 

The  sum  of 

all  wing  and  body  vortex  lattices  may  not  exceed  35. 

Card  : 

Set  8B  -  Vortex  Lattice  Geometry 

Col. 

Variable 

Value 

Description 

1-10 

XBE 

Arbitrary 
(floating  point) 

x  coordinate  of  origin  of 
streamwise  vortex 

11-20 

YBE 

fl 

y  coordinate  of  origin  of  stream- 
wise  vortex 

21-30 

ZBE 

II 

z  coordinate  of  origin  of  stream- 
wise  vortex 

31-40 

MV 

Arbitrary 

Integer 

Number  of  bound  vortices  in  lattice 

2  £  MV  <  60 

41-50 

OPT 

0 

Vortex  lattice  points  to  be  read 
from  card  set  8C 

1 

Vortex  lattice  points  same  as 
preceding.  Omit  card  set  8C. 

2 

Optional  vortex  lattice  control 
panel  coordinates  read  on  card  8D-3. 

51-60 

FLAG 

0 

Normal  case  -  vortex  lattice  panels 
calculated 

1 

Terminal  vortex  of  current  body  vorte: 
lattice  network 

2 

Corner  points  of  control  point  panel 
be  read  on  Cards  8D-2  and  8D-3  (used 

when  arbitrary  control  point  is  desired) 
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Card  Set  8B  -  Vortex  Lattice  Geometry  (cont'd) 


Col.  Variable  Value 


Description 


61-70  SIMPOT  0  Symmetry  option  specified  on  card 

2  enforced  for  this  vortex 

1  Symmetry  option  ignored  for  this 

vortex  lattice  (used  for  inserting 
vortex  lattice  networks  in  vertical 
tails  located  in  x,z  plane) 


Card  Set  8C  -  Vortex  Lattice  Coordinates 


Col. 

Variable 

Value 

Description 

1-10 

B(J) 

Arbitrary 
(floating  point) 

x  coordinate  of  point  J 

11-20 

A(J) 

II 

z  coordinate  of  point  J 

21-30 

C(J) 

II 

Vortex  lattice  strength 

31-40 

D(J) 

II 

Ay  shift  of  point  J 

Card  set  8C  contains  MV  cards  if  OPT  -  0,  and  is  omitted  if  OPT  =  1 
on  card  8B. 

Control  Set  8D  -  Vortex  Lattice  Terminal  Point  &  Control  Point 
Coordinates 


Two  or  three  additional  cards  are  required  to  specify  the  terminal 
point  of  the  streamwise  vortex,  and  the  comer  points  of  the  lattice 
control  point  panel . 

Card  8D-1 


Col. 

Variable 

Value 

Description 

1-10 

B 

Arbitrary 
(floating  point) 

x  coordinate  of  terminal  point 
of  streamwise  vortex 

11-20 

A 

II 

z  coordinate  of  terminal  point 
of  streamwise  vortex 

21-30 

D 

II 

Ay  shift  of  terminal  point  of  stream- 
wise  vortex 

Note:  This  point  also  defines  the  upstream  corner  of  the  control  point 
panel  if  FLAG  ^  2  on  card  8B. 


34 


Card  8D-2  -  Same  as  card  8D-1,  containing  the  coordinates  of  the 
downstream  corner  of  the  control  point  panel  if  FLAG  t  2  on  card  8B. 

If  FLAG  '  2  on  card  8B,  this  card  contains  the  coordinates  of  the 
upstream  corner  of  the  control  point  panel. 

Card  8D-3-  If  FLAG  =  2  on  card  8B,  this  card  contains  the  coordinates 
of  the  downstream  corner  of  the  control  point  panel,  in  the  same  format 
as  card  8D-1.  Omit  this  card  if  FLAG  ^  2  in  card  8B-1. 

Note:  Repeat  card  8B,  and  card  sets  8C  and  8D  NV  times  to  complete 
card  set  8. 

Plot  Input  Cards  -  The  configuration  panel  geometry  is  stored  on 
TAPE  11.  If  PLOT  °  1  on  card  2,  the  plot  overlay  is  called,  and  a 
plot  tape  is  written.  A  sample  panel  geometry  plot  is  shown  on 
Figure  28.  Additional  plot  input  cards  required  are  described  below. 
Omit  these  cards  if  PLOT  =  0  on  card  2. 


Card  9  -  Plot  Parameters 


Col.  Variable  Value 


Description 


1-10  NVU  0 

1-4 


IPRINT  0 

1 

IHIDE  0 

1 

IBUG  0 

1 


No  plot  requested;  return. 

Number  of  view  points  selected.  If 
NVU  positive,  only  source  panels  will 
be  plotted.  If  NVU  negative,  vortex 
panels  will  also  be  plotted. 

Panel  corner  points  are  not  printed. 

Panel  corner  points  printed. 

Eliminate  hidden  surfaces  in  plot 
output . 

All  surfaces  plotted. 

No  debug  printout  from  PLOT 
subroutines. 

Additional  debug  printout  requested. 


Card  Set  10  -  View  Point  Coordinates 


Col. 

Variable 

Value 

Description 

1-10 

VUE(l.NVU) 

Arbitrary 

x  coordinate 

of 

view 

point  NVU 

(floating  point) 

11-20 

VUE(2,NVU) 

ft 

y  coordinate 

of 

view 

point  NVU 

21-30 

VUE(3,NVU) 

If 

z  coordinate 

of 

view 

point  NVU . 
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Card  Set  10  -  View  Point  Coordinates  (cont'd) 


Repeat  Card  Set  10  NVU  times. 

Any  view  point  coordinate  is  set  equal  to  infinity  if  it  is  greater 
than  2^.  (32,768) 

Aerodynamic  Input  Cards 

The  configuration  panel  geometry  is  transferred  to  the  aerodynamic 
section  of  the  program  TAPE  11.  Additional  aerodynamic  input  cards 
required  are  described  below: 

Card  11  -  Case  Identification  Card  -  Card  11  contains  any  desired 


case  Identification 

in  columns  1* 

•80. 

Card  12 

-  Iteration 

Option  Card 

Col. 

Variable 

Value 

Description 

1-10 

NIT 

Arbitrary 

Integer 

Maximum  number  of  iterations 

11-20 

IEPS 

II 

Exponent  of  10  setting  limiting  value 
for  residue  of  iterative  solution  (-3 
or  -4  recommended) 

21-30 

ITYPE 

1 

Gauss-Sledel  iteration  procedure 

2 

Mixed  direct/iterative  solution 
procedure 

Card  13 

-  Configuration  Options 

Col. 

Variable 

Value 

Description 

1-10 

COMPT 

0 

Forces  and  moments  calculated  for 
complete  configuration 

Arbitrary 

Integer 

Forces  and  moments  calculated  on 
components.  Panel  Indices  of  each 
component  follow  on  card  22 

11-20 

SECT 

0 

No  wing  section  forces  and  moments 

1 

Wing  section  forces  and  moments 

calculated.  Wing  section  indices 
follow  on  card  21,  panel  indices  in 
each  section  on  card  22,  and  section 
reference  lengths  on  card  25 
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Card  13  -  Configuration  Options  (cont’d) 


2  Forces  and  moments  calculated  on 

subsections.  The  number  of  sub¬ 
sections  follow  on  card  21,  the 
number  of  panel  groups  on  card  23, 
the  panel  indices  in  each  group 
on  card  24,  and  subsection  reference 
lengths  on  card  25 


Card  14 

-  Reference  Parameters 

Col. 

Variable 

Value 

Description 

1-10 

REFA 

Arbitrary 
(fixed  point) 

Reference  area 

11-20 

REFI 

ft 

Reference  chord  (MAC) 

21-30 

XOO 

II 

Axial  distance  of  leading  edge  of 
MAC  from  origin 

31-40 

X25 

II 

Axial  distance  of  quarter  chord  of 
MAC  from  origin 

Card  15  -  Configuration  Lift  Option 


Col. 

Variable 

Value 

Description 

1-10 

KUT 

0 

Nonlifting  configuration,  no  vortex 
lattice  Kutta  condition  imposed 

1 

Lifting  configuration,  vortex  lattice 
Kutta  condition  imposed 

-1 

Wing  vortex  lattice  extends  through 
body  having  same  strength  as  adjacent 
wing  vortex  lattice 

11-20 

NBV 

Arbitrary 

Integer 

Number  of  body  vortices  (NBV  <  5) 

21-30 

NV  (1) 

II 

Number  of  wing  vortices  associated 
with  body  vortex  1 

31-40 

• 

• 

NV  (2) 

• 

• 

It 

Number  of  wing  vortices  associated 
with  body  vortex  2 

• 

61-70 

NV  (5) 

II 

Number  of  wing  vortices  associated 
with  body  vortex  5 
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Card  16  -  Compressibility  Rule  Option 


Col. 

Variable 

Value 

Description 

1-10 

K0MPR 

1 

Gothert  Rule  1  selected  (See  Text) 

2 

Gothert  Rule  2  selected  (See  Text) 

11-20 

POINTS 

Arbitrary 

Integer 

Number  of  field  points  following  on 
card  set  16A 

21-30 

N0RPAN 

II 

Number  of  normal  velocities  following 
on  card  set  16B 

Card  Set  16A  -  Field  Point  C«  -rdinates 

Col. 

Variable 

Value 

Description 

1-10 

XP 

Arbitrary 
(floating  point) 

Coordinates  of  field  point 

11-20 

YP 

II 

ft  II  II  II 

21-30 

ZP 

II 

If  If  II  II 

(Repeat  card  16A 

POINTS  times.) 

Card  Set  16B  -  Normal  Velocity  Input 

Col. 

Variable 

Value 

Description 

1-10 

NP 

Integer 

Panel  number 

11-20 

NORVEL 

Arbitrary 
(floating  point) 

Normal  velocity  on  panel  NP.  If 

NORVEL  =0,  the  normal  velocity  is  set 
equal  to  the  normal  component  of  the 
onset  veloci ty (i . e .» Equation  A3  is 
used . ) 

Repeat 

card  set  16B  N0RPAN  times. 

Card  Set  17  -  Number  of  Mach  Numbers 

Col. 

Variable 

Value 

Description 

1-10 

NMA 

Arbitrary 

Integer 

Number  of  Mach  numbers  following  on 
tard  set  18 
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Card  Set  18  -  Mach  Number 


Col.  Variable  Value  Description 

1-10  MA  Arbitrary  Mach  number 

(floating  point) 

Card  Set  19  -  Number  of  Angles  of  Attack 
Col.  Variable  Value  Description 

1-10  NAL  Arbitrary  Number  of  angles  of  attack  following 

Integer  on  card  set  20 

Card  Set  20  -  Angle  of  Attack  or  Yaw 

Col.  Variable  Value  Description 

I- 10  ALPHA  Arbitrary  Angle  of  attack  in  degrees 

(floating  point) 

II- 20  BETA  "  Angle  of  yaw  in  degrees 


Card  Set  21  -  Number  of  Sections 

Col.  Variable  Value 

1-10  IS  Arbitrary 

Integer 

Card  Set  22  -  Panel  Indices 
Col .  Variable  Value 

I- 10  IA  Arbitrary 

Integer 

II- 20  IE 

Card  Set  23  -  Number  of  Panels  in 

Col.  Variable  Value 

1-10  IREI  Arbitrary 

Integer 


Description 

Number  of  sections;  omit  if 
SECT  =  0  on  card  13 

Description 

Index  of  initial  panel  in  section 

Index  of  final  panel  in  section 
Subsections 
Description 

Number  of  panels  in  subsections 
Omit  if  SECT  <2  on  card  13 
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Card  Set  24  -  Subsection  Panel  Indices 


Col. 

Variable 

Value 

Description 

1-5 

II  (1) 

Arbitrary 

Integer 

Panel  indices  of  all  panels  in  sub¬ 
section;  omit  if  SECT  <2  on  card  13 

6-10 

II  (2) 

It 

11-15 

II  (3). 

.  .  etc. 

Card 

Set  25  -  Reference  Lengths 

Col. 

Variable 

Value 

Description 

1-10 

DELY 

Arbitrar> 

Integer 

Width  of  section 

11-20 

REEL 

ft 

Reference  length  of  section 

21-30 

XLE 

It 

Moment  reference  point  of  section 

Cards 

21-25  must 

be  repeated  for 

each  angle  of  attack  or  yaw,  if 

section  data  requested. 
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PROGRAM  OUTPUT 


The  standard  output  of  the  program  consists  of  a  list  of  the  input 
cards,  a  table  of  panel  points,  a  table  of  velocities  and  pressure 
coefficients  at  panel  control  points,  and  a  force  and  moment  summary. 
Additional  output  may  be  obtained  by  selecting  appropriate  values  of 
the  integer  PRINT.  A  sample  output  is  given  in  Appendix  IV. 


PRINT  =  1 


PRINT  -  2 


PRINT  -  3 


PRINT  =  4 


Tables  of  panel  corner  points,  centroids,  and 
the  panel  coordinate  transformation  matrix  are 
printed  out.  Individual  panel  forces  and  moments 
are  also  printed  out. 

Tables  of  panel  velocity  components  and  influence 
coefficients  are  printed  out.  This  option  requires 
a  large  line  count  limit. 

The  aerodynamic  influence  coefficient  matrix  is 
printed  out  in  row  order,  together  with  the  right 
side  of  the  matrix  equation,  and  all  solution 
iterations. 

This  option  prints  out  the  successive  solution 
iterations  only. 


PROGRAM  TIME  ESTIMATION 


Estimates  of  the  CPU  time  required  by  the  CDC  6600  computer  to  calculate 
the  aerodynamic  matrix  and  solve  for  one  angle  of  attack  or  yaw  are 
presented  on  Figure  10.  If  the  configuration  is  symmetric  about  the 
x,z  plane  and  the  yaw  angle  is  zero,  the  flow  is  symmetric.  Use  of  this 
fact  by  the  program  reduces  the  running  time.  This  reduction  is  reflected 
in  Figure  10  by  the  fact  that  only  the  number  of  source  and  vortex 
panels  on  one  side  of  the  x,z  plane  are  counted  for  the  symmetric  flow 
case. 

PROGRAM  USAGE 


The  success  of  this  method  of  analysis  depends  to  a  large  extent  on 
the  choice  of  the  number  and  location  of  panels  used  to  represent  the 
configuration.  Certain  features  of  the  program  input  will  be  described 
in  this  section,  together  with  recommendations  on  program  usage. 

Body  Input 

The  body  is  described  by  a  series  of  cross-sections  given  at  selected 
intervals  along  its  length.  The  surface  panels  are  located  between 
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NOTE:  FOR  YAW  CASES  AND  FOR  UNSYMMETRIC  CASES, 
COUNT  ALL  PANELS. 

FOR  SYMMETRIC,  UNYAWED  CASES,  COUNT  PANELS 
ON  ONE  SIDE  OF  SYMMETRY  PLANE  ONLY. 


Figure  10.  CPU  Time  Required  for  CDC  6600. 
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adjacent  sections,  with  the  corner  points  being  defined  by  the  cross- 
section  coordinates.  Unless  the  cross  sections  can  be  described  by 
some  mathematical  formula,  accurate  drawings  are  required  for  each  body 
station.  The  cross  section  can  be  defined  either  in  Cartesion  (y,z) 
or  polar  (r,0)  coordinates  about  the  body  reference  axis.  A  maximum 
of  70  stations  may  be  used  along  the  length  of  the  body,  and  a 
maximum  of  60  points  around  the  half  circumference. 


z 


Body  Section 
MB  =  12 
Origin  at 
(XBE,YBE,ZBE) 


Panel  Corner  Point 


Figure  11.  Body  Cross  Section. 

In  general,  more  panels  are  required  in  regions  with  rapid  changes  of 
cross  sectional  shape,  such  as  around  canopies  or  wing-body  intersections. 

The  program  does  not  require  the  same  number  of  panels  in  each  circumfer¬ 
ential  ring,  and  a  special  input  option  is  provided  to  identify  the 
sections  at  which  the  number  of  panels  is  being  increased  or  decreased. 
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In  addition,  the  panel  corners  may  be  shifted  lengthwise  out  of  the 
plane  of  the  defining  section.  This  option  allows  more  freedom  for 
paneling  complex  wing-body  intersections  and  fairings. 

Wing  Input 

The  wing  is  described  by  a  series  of  airfoil  sections  given  at  selected 
intervals  along  the  span.  The  surface  panels  are  located  between 
adjacent  sections,  with  the  corner  points  being  defined  by  the  section 
coordinates.  The  section  coordinates  may  be  given  in  percent  chord  or 
directly  in  terms  of  the  reference  coordinate  system.  The  panels  in 
each  wing  section  are  generally  numbered  sequentially  from  the  trailing 
edge  on  the  lower  surface  around  the  leading  edge  to  the  trailing  edge 
on  the  upper  surface.  The  same  number  of  points  at  approximately  the 
same  percent  chord  locations  must  be  used  to  define  the  wing  upper  and 
lower  surfaces,  since  the  vortex  lattice  panels  on  the  mean  camber 
surface  are  defined  by  averaging  the  upper  and  lower  surface  points. 

A  maximum  of  60  points  may  be  used  to  define  each  section,  and  a  maximum 
of  40  sections  may  be  defined  on  the  half-wing. 


Wing  Section 
MW  =  21 


Figure  12.  Wing  Section. 
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The  program  does  not  require  the  same  number  of  points  on  adjacent  wing 
sections,  and  a  special  input  option  is  provided  to  identify  the  sections 
at  which  the  number  of  panels  is  being  increased  or  decreased.  In 
addition,  the  panel  corners  may  be  shifted  spanwise  out  of  the  plane  of 
the  defining  section  to  aid  in  the  paneling  of  wing  tips  or  complex 
wing-body  intersections.  In  general,  more  panels  are  required  in  regions 
of  rapid  curvature,  for  example,  the  leading  edge  region. 

Wing  tip  paneling  may  be  omitted  for  wings  having  a  maximum  thickness 
of  less  than  5%.  However,  wing  tip  panels  can  be  included  as  special 
body  panels,  or  read  in  individually  using  the  single  panel  input 
option. 

Vortex  Lattice  Panels 


Vortex  lattice  panels  may  be  placed  inside  either  the  wing  or  body, 
or  omitted  entirely.  For  configurations  with  a  wing,  the  vortex 
lattice  panels  are  automatically  located  on  the  mean  camber  surface 
of  the  wing.  The  vortices  extend  a  finite  distance  behind  the  wing 
in  a  plane  passing  through  the  trailing  edge  and  bisecting  the  trailing 
edge  angle.  The  vortices  should  be  allowed  to  extend  at  least  ten  chord 
lengths  behind  the  wing  to  give  a  reasonable  approximation  of  the  wake. 

A  control  panel  is  associated  with  each  vortex  lattice  and  sized  such 
that  the  panel  control  point  is  located  one  percent  of  the  local  chord 
behind  the  trailing  edge  of  the  wing. 

The  relative  strengths  y^  of  the  individual  bound  vortices  making  up 

the  vortex  lattices  must  be  specified  in  advance.  It  is  recommended 
that  these  strengths  be  proportional  to  the  airfoil  thickness  at  the 
chordwise  location  of  the  bound  vortex.  The  accuracy  of  the  final 
solution  depends  to  some  extent  on  the  vortex  distribution  selected, 
so  some  adjustment  to  the  bound  vortex  strengths  may  be  necessary  if 
a  poor  initial  choice  has  been  made. 

The  program  requires  that  the  y  array  be  read  in  order  of  increasing 
chordwise  station.  Since  the  wing  numbering  system  starts  at  the 
trailing  edge,  the  first  (M+l)/2  points  are  set  equal  to  zero,  and  the 
desired  y  array  associated  with  the  remainder. 

The  structure  of  a  vortex  lattice  is  illustrated  in  the  following 
sketch. 


45 


Leading 

Edge  Trailing 


Control  Point 


Figure  13.  Wing  Vortex  Lattice. 

It  can  be  seen  that  the  trailing  vortices  in  the  wake  are  made  up 
of  the  sum  of  the  streamwise  legs  of  the  individual  bound  vortices. 

Vortex  lattices  may  also  be  added  inside  the  body  to  provide  a 
mechanism  for  generating  body  lift.  For  wing-body  combinations,  a 
special  option  is  provided  to  give  the  body  vortex  lattice  the  same 
strength  as  the  adjoining  wing  vortex  lattice,  so  that  the  inboard 
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trailing  vortex  from  the  wing  vortex  lattice  will  be  exactly 
cancelled  by  the  outboard  trailing  vortex  from  the  body  vortex  lattice. 


Figure  14.  Vortex  Lattice  Inside  Body, 


Body  vortices  are  also  used  to  generate  the  circulation  about  vertical 
tiils  located  in  the  plane  of  symmetry  for  yawed  configurations.  This 
technique  must  be  employed  since  the  wing  vortex  lattices  defined  by 
the  program  will  automatically  cancel  in  the  plane  of  symmetry.  For 
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unyawed  configurations,  however,  no  wing  or  body  vortices  may  be  used 
in  vertical  tails  located  in  the  plane  of  symmetry.  The  use  of  body 
vortices  in  a  vertical  tail  is  illustrated  in  Figure  15. 


Figure  15.  Vortices  in  Vertical  Tail. 


For  configurations  having  wings,  horizontal  tails,  and  vertical 
tails,  and  employing  body  vortices  to  carry  the  lift  generated  by  the 
wing  and  horizontal  tail  through  the  body,  the  body  vortices  used  to 
provide  lift  on  the  vertical  tail  must  be  input  after  the  body  vortices 
associated  with  the  wing  or  horizontal  tail. 
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A  GENERAL  GUIDE  TO  PANELING  WING-BODY-NACELLE  CONFIGURATIONS 


An  input  option  is  available  to  take  advantage  of  symmetry.  With  the 
x  z  plane  in  the  plane  of  symmftry,  only  the  half  of  the  configuration 
corresponding  to  positive  y  is  required  as  input.  This  rule  applies 
whether  or  not  the  free  stream  is  in  the  plane  of  symmetry  (i.e.,  with 
or  without  yaw) . 

When  paneling  a  particular  configuration,  attention  should  be  paid  to 
the  principal  regions  of  interest,  since  this  allows  for  optimum  use 
of  the  total  number  of  panels.  A  denser  panel  distribution  should  be 
used  in  the  region  of  interest  and  a  lesser  number  of  panels  in  regions 
removed  from  this  region.  A  sparse  panel  distribution  has  a  substantial 
local  effect.  When  choosing  the  number  of  panels  to  be  used  in  each 
region,  it  is  important  to  realize  that  the  pressure  calculated  at  the 
centroid  of  each  panel  is  assumed  to  exist  over  the  entire  panel.  There¬ 
fore,  more  panels  are  required  in  regions  of  large  pressure  gradients. 

The  three-dimensional  potential  flow  program  usually  accepts  panel 
corner  points  as  inputs.  An  exception  to  this  is  the  case  where  a 
particular  section  is  defined  as  a  circular  arc,  which  may  then  be 
specified  along  with  the  number  of  equispaced  panel  corner  points  along 
the  circle.  The  program  accepts  the  defined  panels  in  networks  -  each 
network  defining  a  specific  region  of  the  configuration.  The  networks 
when  assembled  together  define  the  complete  configuration. 

In  setting  up  the  complete  paneling  scheme  for  a  wing-body-nacelle 
configuration,  it  is  usually  best  to  define  the  nacelle,  wing,  and 
body  as  separate  sets  of  networks.  The  regions  where  two  bodies  inter¬ 
sect  are  then  treated  separately,  and  the  particular  networks  require 
redefining.  Any  fairings  that  exist  are  treated  as  part  of  the 
particular  configuration  and  not  as  a  separate  configuration.  Each 
fairing  is  usually  specified  by  its  own  network  (or  set  of  networks) 
of  panels. 

The  body,  or  most  likely  half-body,  is  usually  paneled  by  specifying 
the  buttock  line  (y)  and  waterline  (z)  coordinates  at  each  body  station 
(x) ,  such  that  at  each  body  station  the  panel  points  are  equispaced. 

(It  is  not  necessary  to  have  the  panels  equispaced,  but  it  is  often 
more  convenient.) 

The  nacelle  is  paneled  in  a  manner  similar  to  the  body. 

The  wing  is  usually  paneled  by  specifying  the  x  and  y  coordinates  at 
each  buttock  line  (y) .  When  defining  the  panels,  it  must  be  remembered 
that  more  panels  are  required  at  the  leading  edge,  where  the  pressure 
gradients  are  large.  The  panel  size  should  not  vary  by  more  than  50% 
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(larger  or  smaller)  from  any  adjacent  panel.  It  is  common  practice 
to  have  no  panel  larger  than  5  percent  of  the  wing  chord.  Constant 
spacing  is  the  optimum  scheme  for  spanwise  paneling. 

The  paneling  at  the  intersection  of  two  segments  requires  special 
care.  For  the  case  of  wing-body  intersection,  the  body  panels  above 
and  below  the  wing  must  be  adjusted  to  account  for  the  area  eliminated 
by  the  intersection.  At  each  body  station  (x) ,  the  y  and  z  coordinates 
are  adjusted  to  give  approximate  equispacing  above  and  below  the  inter¬ 
section.  For  the  wing  segment,  the  existing  paneling  can  be  maintained. 
All  intersections  are  treated  in  this  manner. 

The  wing  requires  a  system  of  multihorseshoe  vortices,  placed  inside 
the  wing  with  the  trailing  vortices  emanating  from  the  trailing  edge, 
to  produce  lift.  The  number  of  chordwise  locations  of  the  internal 
bound  vortices  are  chosen  to  minimize  large  local  disturbances  at 
the  wing  surface.  The  vortices  are  placed  along  the  camberline, 
equidistant  from  the  nearest  panel  corner  points  in  the  chordwise 
direction. 

Additional  segments  may  be  added  to  the  configuration,  such  as  a 
vertical  tail  and  a  horizontal  tail.  The  intersections  of  the  various 
segments  are  treated  as  above;  and  both  horizontal  and  vertical  tails 
will  require  internal  vortex  lattices  to  produce  lift.  No  vortex 
lattices  are  placed  in  vertical  tails  located  in  the  plane  of  symmetry 
unless  the  configuration  is  yawed. 

It  should  be  noted  that,  except  for  the  simplest  configurations,  the 
preparation  of  the  input  to  the  three-dimensional  potential  flow 
program  is  cumbersome  and  usually  almost  always  requires  the  aid  of 
auxiliary  geometry  manipulation  computer  programs.  The  use  of  the 
three-dimensional  plotting  program  is  also  almost  essential  to  check 
the  panel  input  data. 

Save  Tape  Option 

Since  a  major  portion  of  the  computer  time  is  taken  up  by  the  calculation 
of  the  aerodynamic  matrices,  provision  is  made  for  saving  these  matrices 
on  magnetic  tape  for  subsequent  runs  on  the  same  configuration.  A  new 
tape  must  be  generated  for  each  Mach  number,  however.  The  program 
stores  the  aerodynamic  matrix  on  auxiliary  disc  file  TAPE  10,  the 
geometric  data  on  auxiliary  disc  file  TAPE  11,  and  the  velocity 
component  matrices  on  auxiliary  disc  file  TAPE  12.  If  the  save  tape 
option  is  selected,  ISAVE  =  -1,  a  magnetic  tape  must  be  designated 
to  replace  the  disc  file  TAPE  11.  The  contents  of  TAPE  10  and  TAPE 
12  are  also  transformed  to  this  tape  during  the  run. 

On  subsequent  runs,  the  contents  of  the  tape  must  be  transferred  back 
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to  the  three  disc  files,  and  the  program  rerun  with  ISAVE  =  1.  Only 
the  first  two  cards  from  the  geometry  input  and  the  aerodynamic  input 
(Cards  11-25)  are  required  if  this  option  is  selected. 
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COMPARISON  WITH  EXPERIMENT 


The  purpose  of  this  section  is  to  aid  the  reader  in  the  evaluation  of 
the  computer  program.  Discussion  of  comparisons  between  theory  and 
experiment  will  be  limited  to  two  specific  configurations,  although 
many  configurations  have  been  investigated  at  various  times  using 
computer  program  WBAERO. 

BO  105  HELICOPTER  CONFIGURATION 


The  geometry  analyzed  is  shown  in  its  paneled  configuration  in  Figure 
16.  The  number  of  panels  on  one  side  of  the  plane  of  symmetry  is  256. 

The  results  for  calculations  with  separation  modeling  were  given  eariler 
(see  Figures  7  and  8)  in  the  explanation  of  this  modeling.  Calculations 
with  no  separation  modeling  are  compared  with  experimental  data  in 
Figures  17  through  20  for  the  cases  of  angle  of  attack  ot  =  0°,  and  angle 
of  yaw  6  -  0  ,  as  well  as  a  -  0  t  8  =  *  10.  Calculations  have  also 
been  made  by  Gillespie  (8)  for  the  same  configuration  using  the  Douglas 
Neumann  program,  and  these  have  been  included  for  comparison  purposes. 

In  Figure  17,  calculated  and  measured  pressure  coefficients  are  shown 
as  a  function  of  axial  distance  for  the  BO  105  fuselage  top  centerline 
with  a  =  0  ,8  =  0  .  Since  the  Douglas  Neumann  and  WBAERO  programs 
are  fundamentally  identical  for  the  nonlifting  case,  it  would  be  expected 
that  the  calculated  results  should  be  identical.  On  closer  scrutiny, 
however,  it  is  realized  that  the  programs  actually  differ  in  two 
respects.  First,  WBAERO  used  the  exact  expression  for  the  velocity 
perturbation  due  to  a  source  for  a  greater  distance  away  from  the  source 
panel  than  does  the  Douglas  Neumann  program,’  and  second,  the  iterative 
techniques  used  by  each  program  in  inverting  the  influence  coefficient 
matrix  are  different.  While  the  latter  difference  does  not  necessarily 
affect  the  accuracy  of  solution,  the  employment  of  different  expressions 
for  the  velocity  perturbation  will  most  certainly  affect  accuracy,  and  it 
is  reasonable  to  expect  WBAERO  to  be  slightly  more  accurate  than  the 
Douglas  Neumann  program.  Such  appears  to  be  the  case  for  the  data  of 
Figure  17,  at  least  until  the  region  of  flow  separation  is  approached. 
Comparisons  are  shown  for  the  same  angles  of  attack  and  yaw  in  Figure 
18  for  pressure  coefficients  along  Waterline  10  (see  Figure  16).  As 
in  Figure  17,  both  methods  are  in  good  agreement  with  experiment.  A 
further  possible  reason  for  the  slight  discrepancy  between  the  two 
calculation  procedures  results  from  the  panel  distribution  used  to 
represent  the  geometry.  Identical  panel  distributions  were  used  in 
each  program  up  to  axial  station  32  {  from  that  point  on  slightly 
different  panel  arrangements  were  used  to  represent  the  body  closure. 
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Additional  calculations  were  made  for  comparison  with  measurements 
at  a  =  0  t  0  =10  .  As  shown  in  Figure  19,  agreement  between  theory 
and  experiment  for  pressures  along  the  top  centerline  is  quite  good. 
Comparisons  between  theory  and  experiment  for  both  windward  and  leeward 
sides  of  the  BO  105  fuselage  at  Waterline  6  show  good  agreement  as  seen 
on  Figure  20. 

HEAVY  LIFT  HELICOPTER  (HLH)  CONFIGURATION 

This  configuration  is  shown  as  paneled  in  Figure  21.  The  number  of 
panels  on  one  side  of  the  plane  of  symmetry  is  665.  In  the  analysis 
of  this  configuration  the  boundary  condition  discussed  in  the  section 
on  separation  modeling  was  used  in  order  that  large  suction  pressure 
peaks  could  be  avoided  at  sh_.p  corners.  In  regions  immediately  behind 

these  corners  the  condition  V  •  n  is  everywhere  zero  has  been  relaxed  such 

that  V  •  n  =  V  •  n.  Separated  flow  regions  on  the  wings  as  on  the 

00 

nacelle  struts  have  not  been  modelled  at  this  time.  Calculations  have 
been  performed  for  zero  yaw  at  two  angles  of  attack,  ct  =  0,  -8  ,  The 
results  of  this  analysis  are  compared  with  experimental  data  obtained 
from  Reference  9  at  various  locations  on  a  1/12  scale  HLH  configuration 
as  shown  in  Figures  22  through  27.  Theoretical  and  measured  pressure 
coefficients  along  the  front  pylon  top  centerline  are  shown  in  Figure 
22.  Agreement  between  theory  and  experiment  is  quite  good  in  the  region 
ahead  of  axial  station  10.  At  approximately  that  location  the  wind 
tunnel  model  has  an  attachment  for  the  forward  hub  which  was  not  modelled 
in  the  potential  flow  calculations.  Consequently,  the  calculated 
pressures  are  not  expected  to  be  in  close  agreement  with  experiment  in 
the  local  region  behind  the  hub  attachment.  Calculations  along  the 
front  pylon  bottom  centerline  (Figure  23)  give  a  good  indication  of 
the  effect  of  sharp  corner  modelling  on  the  comparisons  with  experiments. 
In  this  case  the  calculation  by  Gillespie  using  the  XYZ  program 
(Reference  7)  shows  a  large  suction  peak  around  the  shaip  corner  of  the 
fuselage  leading  to  the  hoist  operator's  observation  window.  The 
calculations  using  WBAER0  while  showing  some  overshoot  in  pressure  are 
in  much  better  agreement  with  experiment,  while  at  the  same  time  taking 
less  computer  time  because  of  the  more  accurate  modelling  of  the  real 
flow.  Figures  24  and  25  give  the  results  of  comparisons  between  theory 
and  experiment  for  the  wing  upper  and  lower  surface  pressures.  Although 
separation  effects  were  not  modelled  except  behind  sharp  corners,  WBAER0 
is  generally  in  good  agreement  with  experiment.  Further  comparisons  have 
been  made  for  the  nacelle  at  the  maximum  spanwise  location.  Again, 
the  effect  of  a  sharp  corner  has  been  accounted  for,  with  excellent 
results  as  shown  on  Figure  26.  Measured  and  calculated  pressure 
coefficients  have  also  been  compared  along  waterlines  for  the  aft 
pylon.  In  Figure  27,  the  comparison  is  for  a  waterline  just  above  the 
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nacelle  strut.  One  particularly  interesting  result  of  this  comparison 
is  the  suction  peak  resulting  from  interference  by  the  nacelle  strut 
which  is  predicted  by  theory.  Because  of  the  choice  of  pressure  tap 
locations,  no  indication  of  this  peak  can  be  obtained  from  data 
alone.  If  load  calculations  are  made  from  such  measurements,  erroneous 
conclusions  can  be  drawn. 

SAMPLE  CASE 


A  simple  wing-body-vertical-tail  airplane  configuration  (Figure  28) 
has  been  analyzed  as  a  sample  case  in  order  to  demonstrate  most  of 
the  program  features.  Because  configurations  having  large  vertical 
tails  or  pylons  are  expected  to  experience  side  forces  under  yawing 
conditions,  it  was  necessary  to  modify  the  computer  program  in  order 
that  vortex  networks  could  be  employed  to  provide  the  correct  circulation 
for  side  force  calculations.  Calculations  are  shown  in  Figure  29  for 
the  pressure  distribution  along  a  waterline  on  the  vertical  tail.  Two 
calculations  are  given,  for  the  conditions  8  =  0  and  8  =  10  .  It  is 
interesting  to  note  that  for  8  =  10  the  vertical  tail  carries  a  nega¬ 
tive  load  on  the  aft  portion.  Slender  wing  theory  for  low  aspect 
ratio  wings  suggests  that  there  should  be  no  load  over  the  aft  part 
of  the  wing;  however,  interference  with  the  blunt  based  fuselage 
appears  to  lead  to  the  negative  load. 
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Figure  17.  Pressure  Distribution  for  BO  105  along  Fuselage 
Top  Centerline  a  =  0°  g  =  0° • 
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Figure  18.  Pressure  Distribution  for  BO  105  along  Fuselage 
Waterline  10  a  =  0°.  8=0°. 
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Figure  19.  Pressure  Distribution  for  BO  105  along  Fuselage 
Top  Centerline  a  =  0°,  g  =  10°. 
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Pylon  Top  Centerline. 
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Pylon  Bottom  Centerline. 
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Figure  24.  Pressure  Distribution  for  HLH  Wing  Upper  Surface 
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Figure  25.  Pressure  Distribution  for  HLH  Wing  Lower  Surface 
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Figure  26.  Pressure  Distribution  for  HLH  Nacelle  -  Maximum  Span. 
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Figure  27.  Pressure  Distribution  for  HLH  Aft  Pylon  Above 
Nacelle  Z  =  11.3. 
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Combination. 


Figure  29.  Pressure  Distribution  on  a  Low  Aspect  Ratio 
Vertical  Tail. 


CONCLUSIONS 


As  a  result  of  the  studies  described  in  this  report,  the  following 
conclusions  have  been  drawn: 

1.  The  three-dimensional  potential  flow  computer  program  provides 
the  user  with  a  practical  and  accurate  method  for  the  calculation 
of  pressures  and  aerodynamic  forces  for  arbitrary  shaped  lifting 
configurations  including  the  effect  of  yaw. 

2.  The  separation  modelling  proposed  in  this  report  shows  encourag¬ 
ing  agreement  with  experiment  while  reducing  the  computer  time 
requirements  for  a  given  configuration. 
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APPENDIX  I 


PANEL  GEOMETRY  CALCULATIONS 


The  analytical  procedure  presented  here  follows  closely  the  method 
first  developed  in  Reference  5.  A  quadrilateral  surface  element  is 
described  by  four  corner  points,  not  necessarily  lying  in  the  same 
plane,  as  shown  in  the  sketch.  The  quadrilateral  element  is  approxi¬ 
mated  by  a  planar  panel  as  follows: 


The  coordinates  in  the  reference  coordinate  system  are  identified  by 

-►  ■> 

their  subscripts.  The  components  of  the  diagonal  vectors  T^  and  T2  are 

*1*  -  *3  -  Xj  tly  -  y3  -  ?1  tu-  s3  -  ^ 

C2x  '  *4  -  *2  C2y  '  *4  -  *2  l2z’  z4  ‘  z2 

We  may  now  obtain  a  vector  N  (and  its  components)  by  taking  the  cross 
product  of  the  diagonal  vectors. 


-y  •* 

N  *  t2  x 


N  - 

X 

t2ytlz 

t,  t_ 
ly  2z 

N  - 

y 

tlxt2z 

t2xtlz 

Nz  - 

t2xtly 

tlxt2y 
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The  unit  normal  vector,  n,  to  the  plane  of  the  element  is  taken  as  N 

divided  by  its  own  length  |n|  (direction  cosines  of  outward  unit 
normal) . 


where 


N 

N 

-JL 

N 

N 

_ z 

N 


it 


+  N 


N2]1/2 

z 


The  plane  of  the  element  is  now  completely  determined  if  a  point  in  this 
plane  is  specified.  This  point  is  taken  as  the  point  whose  coordinates 

x,  y,  z  are  the  average  of  the  coordinates  of  the  four  input  points. 


x 

y 

z 


4  [X1 

+  X2 

+  X3 

+ 

x4] 

£  [yi 

+  y2 

+  y3 

+ 

£  [zi 

+  Z2 

+  Z3 

+ 

Z4] 

Now  the  input  points  will  be  projected  into  the  plane  of  the  element 
along  the  normal  vector.  The  resulting  points  are  the  corner  points  of 
the  quadrilateral  element.  The  input  points  are  equidistant  from  the 
plane,  and  this  distance  is 


d 


I  nx(x  -  xx)  +  ny(y 


yx>  +  nz(z  -  zx)  | 
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The  coordinates  of  the  corner  points  In  the  reference  coordinate  system 
are  given  by 


» 


,  . .  k+1  , 

+  ('1)  V 

yk  +  (-l)k+1nyd 
,  ,  v  k+1 

zk  +  (_1)  nzQ 


k 


1,  2,  3,  A 


Now  the  element  coordinate  system  must  be  constructed.  This  requires 
the  components  of  three  mutually  perpendicular  unit  vectors,  one  of 
which  points  along  each  of  the  coordinate  axis  of  the  system,  and  also 
the  coordinates  of  the  origin  of  the  coordinate  system.  All  these 
quantities  must  be  given  in  terms  of  the  reference  coordinate  system. 

The  unit  normal  vector  is  taken  as  one  of  the  unit  vectors,  so  two 
perpendicular  unit  vectors  in  the  plane  of  the  element  are  needed.  Denote 

these  unit  vectors  t^  and  The  vector  t^  is  taken  as  T^  divided  by 

its  own  length  T^,  i.e., 


t 


lx 


t 


ly 


T 


1 


T 


1 


t 


lz 


T 


1 


where 


[t: 


lx 


+  T 


iy 


tL]1/2 
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The  vector  is  defined  by  -  n  *  t^,  so  that  its  components  are 


2x 

=  n  t, 
y  lz 

-  n  t. 
z  ly 

2y 

=  n  t, 
z  lx 

-  n  t. 
x  lz 

2z 

=  n  t. 
x  ly 

-  n  t. 
y  lx 

The  vector  t^  is  the  unit  vector  parallel  to  the  x  or  E,  axis  of  the 
element  coordinate  system,  while  tj  is  parallel  to  the  y  or  ri  axis,  and 
n  is  parallel  to  the  z  of  ?  axis  of  this  coordinate  system. 


To  transform  the  coordinate  of  points  and  the  components  of  vector 
between  the  reference  coordinate  system  and  the  element  coordinate 
system,  a  transformation  matrix  is  required.  The  elements  of  this  matrix 

are  the  three  components  of  the  three  basic  unit  vectors  t^,  tL,  and  n. 


T 


lx 

tiy 

fclz 

2x 

S 

t2z 

t 

n 

n 

X 

y 

z 

In  the  computer  program,  the  elements  of  this  matrix  are  referred  to 
as  follows: 


all  = 

fclz 

a12 

ciy 

a13 

=  tlz 

a21  = 

l2x 

a22 

fc2y 

a23 

t2z 

a31  = 

n 

X 

a  32 

=  n 

y 

a33 

=  n 

z 

The  corner  points  are  now  transformed  into  the  element  coordinate 
system  based  on  the  average  point  as  origin.  These  points  have  coor- 

*  i  i 


dinates  xk,  yk, 


zk  in  the  reference  coordinate  system. 


Their  coordinates 


in  the  element  coordinate  system  with  this  origin  are  denoted  by  £k,  rik, 
0.  Because  they  lie  in  the  plane  of  the  element,  they  have  a  zero  z  or  ? 
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coordinate  in  the  element  coordinate  system.  Also,  because  the 

vector  t^,  which  defines  the  x  or  5  axis  of  the  element  coordinate 

system,  is  a  multiple  of  the  "diagonal"  vector  from  point  1  to  3,  the 
coordinate  and  the  coordinate  are  equal.  In  the  (£,ri)  coordirate 
system,  the  corner  points  of  the  element  are: 


e 


k 


\  * 


(i  '  Xi>  +  V*  'V  +  hz 

<*  -  V  +  V*  -yl?  +  C2z 


(z  -V 

u  -z'k) 


These  corner  points  are  taken  as  the  comers  of  a  plane  quadrilateral  as 
illustrated  in  the  following  sketch. 


The  origin  of  the  element  coordinate  system  is  now  transferred  to  the 
centroid  of  the  area  of  the  quadrilateral.  With  the  average  point  as 
origin  the  coordinates  of  the  centroid  in  the  element  coordinate  system 
are: 


-  T  n4  u«(ni  -  V  +  «2<1«  -  V> 
"o  *  y11! 
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These  are  subtracted  from  the  coordinates  of  the  comer  points  in  the 
element  coordinate  system  based  on  the  average  point  as  origin  to  obtain 
the  coordinates  of  the  corner  points  in  the  element  coordinate  system 
based  on  the  centroid  as  origin.  Accordingly,  these  latter  coordinates 
are 


k  =  1,  2,  3,  4 

\  "  n0 


Since  the  centroid  is  to  be  used  as  the  control  point  of  the  element, 
its  coordinates  in  the  reference  coordinate  system  are  required.  These 
coordinates  are 


x0  = 

*  +  Vo  + 

Vo 

y0  = 

y  +  cly50  + 

Vo 

z0  " 

5  +  tlz50  + 

Vo 

Finally,  the  area  of  the  quadrilateral  is 

A  *=  j  (S3  -  5j)(n2  -  n4) 
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APPENDIX  II 


SUBROUTINE  DESCRIPTIONS 
(arranged  in  alphabetical  order) 


This  appendix  contains  a  brief  outline  of  the  purpose,  method,  and 
use  or  each  subroutine.  The  principal  constants  and  variables  in 
each  are  listed  in  the  order  of  their  first  appearance,  and  identified 
as  input  or  output  data. 
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Subroutine  ANALOG 


Purpose:  To  transform  the  panel  control  point  and  corner  points 

from  the  real  body  to  the  analog  body,  as  required  by 
Gothert's  compressibility  rule. 

Method:  The  panel  control  point  is  transformed  simply  by  multi¬ 

plying  the  y  and  z  coordinates  by  B  =  -J l-M^.  The  trans¬ 
formation  matrix  is  then  corrected  for  compressiblity,  and 
new  panel  corner  points  determined  in  the  panel  reference 
system.  Finally,  the  area  and  maximum  diagonal  of  the 
panel  are  calculated. 

Use:  Call  ANALOG  (BETA,  YC,  ZC,  A,  XEI,  YEI ,  F,  T) 


2 

Prandtl-Glauert  factor  =  1-M 

y  coordinate  of  panel  control  point 

ZC  z  coordinate  of  panel  control  point 

A  Transformation  matrix 


XEI  Transformed  x  coordinate  of  panel  corner  point  in 

panel  coordinate  system 

YEI  Transformed  y  coordinate  of  panel  corner  point  in 

panel  coordinate  system 

F  Panel  area 
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Subroutine  ANALOG(cont ' d ) 


Output : (cont'd) 

T  Length  of  maximum  diagonal  of  panel 

Subroutine  Called:  None 

Error  Returns:  None 
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Subroutine  DRAW3D 


Purpose: 
Method : 


Use: 


To  plot  the  object  as  viewed  from  the  viewpoints. 

DRAW3D  begins  by  projecting  the  object  onto  the  display 
plane  (the  "V,  W"  plane)  perpendicular  to  the  vector  from 
the  origin  to  the  viewpoint.  It  also  generates  the  panels 
needed  for  symmetry  if  ISYM  -t.s  zero.  The  V,  W  data  are 
then  scaled  so  that  they  will  be  plotted  in  an  area  6.5 
inches  wide  by  9.0  inches  high. 

The  program  then  creates  a  new  list  of  the  panels  to  be 
plotted  in  the  array  LL.  If  the  IHIDE  parameter  is  zero, 
the  new  list  will  only  have  the  panels  that  face  the 
viewpoints,  and  in  addition  the  subroutine  OVRLAP  will 
be  called  to  eliminate  any  panels  that  are  partially  or 
entirely  blocked  by  some  other  panel. 

Finally,  the  program  loops  through  the  new  list  (LL)  and 
calls  the  subroutine  PLTPAN  to  plot  the  panels. 

Call  DRAW3D  (N,  M,  VX,  VY,  VZ,  ISYM,  IHIDE,  IBUG) 

N  The  number  of  vertices 

M  The  number  of  panels 

VX,VY,  The  coordinates  of  the  viewpoint 

VZ 

ISYM  If  zero,  the  object  i6  symmetric  about  the 
XZ  plane 

IHIDE  If  zero,  the  program  is  to  eliminate  the  hidden 
surfaces 

IBUG  If  nonzero,  the  coordinate  of  the  vertices  in  the 
display  plane  will  be  printed.  If  IBUG  =  1,  the 
list  contained  in  the  array  LL  will  be  printed 
along  with  the  corresponding  coordinates  in  the 
display  plane 

Input : 

In  addition  to  the  parameters  in  the  calling  sequence,  the 
subroutine  DRAW3D  uses  the  data  in  the  common  arrays  V,  W, 
DIST,  and  L.  These  arrays  are  set  up  in  the  program  WBPLOT. 
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Subroutine  DRAW3D(cont ' d ) 


Output : 

The  primary  output  of  DRAW3D  is  the  data  in  the  common 
arrays  V,  W,  DIST,  and  LL,  which  arc  used  by  the  subroutine 
OVRLAP  and  PLTPAN.  In  addition  to  this  there  are  the 
optional  IBUG  printouts,  and  label  on  the  plot. 

Subroutines  Used: 

NUMBER  Used  to  draw  numbers  on  the  plot 
OVRLAP  Used  to  eliminate  blocked  panels 
PLOT  Used  to  control  the  plotting  origin 
PLTPAN  Used  to  draw  the  panels 

rvMur)L  Used  to  draw  symbols  on  the  plot 

Limitations:  See  WBPLOT 
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Subroutine  ELEMEN 


Purpose : 

Method : 
Use : 


To  calculate  the  panel  control  point,  corner  points,  and 
transformation  matrix. 


The  method  is  described  in  Appendix  I  of  this  report. 
Call  ELEMEN  (  I,  X,  Y,  Z,  KUTA,  XC,  YC,  ZC,  XE,  YE,  A) 
Input : 

I  Panel  number  (not  used) 


KUTA 

Output : 

XC 

YC 

ZC 


XE 


YE 


Subroutine  Called: 
Error  Returns: 


x  coordinate  of  panel  in  reference  coordinate 
system 

y  coordinate  of  panel  in  reference  coordinate 
system 

z  coordinate  of  panel  in  reference  coordinate 
system 

Not  used 


x  coordinate  of  panel  control  point  (centroid)  in 
reference  coordinate  system 

y  coordinate  of  panel  control  point  (centroid)  in 
reference  coordinate  system 

z  coordinate  of  panel  control  point  (centroid)  in 
reference  coordinate  system 

x  coordinate  of  panel  control  point  in  panel 
coordinate  system 

y  coordinate  of  panel  corner  point  in  panel 
coordinate  system 

Transformation  matrix  (to  transform  from  reference 
coordinate  system  to  panel  coordinate  system) 


None 

None 
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Subroutine  FORMOM 


Purpose: 

Method : 


Use: 


To  calculate  the  force  and  moment  coefficients  on  wing 
and  body. 

The  panel  area,  control  point  coordinates,  and  direction 
cosines  are  read  from  TAPE  13  for  each  panel  in  sequence. 
The  normal  force,  lateral  force,  axial  force,  and  pitching 
moments  of  the  panel  about  the  origin  of  coordinates  are 
calculated,  and  summed.  The  total  force  and  moment 
coefficients  acting  on  the  configuration  are  then  obtained 
by  dividing  these  sums  by  the  reference  area.  The  lift, 
side  force,  and  drag  are  obtained  by  resolving  the  normal 
force,  lateral  force,  and  axial  force  coefficients  into 
the  wind  axis  system. 

If  wing  section  data  is  required,  the  program  calculates 
the  spanwise  lift  and  drag  distribution  on  the  wing.  The 
forces  and  moment  acting  on  each  panel  in  a  given  column 
are  summed,  and  the  section  force  and  moment  coefficients 
are  calculated  with  reference  to  the  area  of  the  column 
of  panels. 

The  output  is  summarized  in  tables  giving  the  input 
geometrical  data,  pressure  coefficients,  panel  forces  and 
moment.  The  integrated  force  and  moment  data  on  the  wing 
sections  and  the  complete  configuration  are  also  tabulated. 

Call  FORMOM,  (ALPHA,  BETA,  MA,  NPAN,  SECT,  COMPT,  REFA, 
REFL,  X00,  X25) 

Input : 

NPAN  Number  of  panels  on  the  configuration 

MA  Mach  number 

ALPHA  Angle  of  attack  (degrees) 

BETA  Angle  of  yaw  (degrees) 

SECT  Wing  section  data  parameter 

COMPT  Component  indicator 

REFA  Reference  area 

REFL  Reference  length 
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Subroutine  F0RM0M(cont 'd) 

Input: 

(cont'd) 

XOO 

Distance  of  leading  edge  of  MAC 

from  origin 

X25 

Distance  of  quarter  chord  of  MAC 

from  origin 

CP 

Pressure  coefficient 

F 

Panel  area 

F1.F2, 

Direction  cosines  of  normal 

F3 

XC,YC, 

Panel  centroid  coordinates 

ZC 

DELY 

Width  of  column  of  panels 

XLE 

Distance  from  column  leading  edge  to  origin 

Output 

DCZ 

Panel  normal  force 

DCY 

Panel  lateral  force 

DCX 

Panel  axial  force 

DCMX 

Panel  moment  about  x  axis 

DCMY 

Panel  moment  about  y  axis 

DCMZ 

Panel  moment  about  z  axis 

CZ 

Normal  force  coefficient 

CY 

Lateral  force  coefficient 

CX 

Axial  force  coefficient 

CL 

Lift  coefficient 

CS 

Side  force  coefficient 

CD 

Drag  coefficient 

CMX 

Pitching  moment  about  x  axis 
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Subroutine  F0RM0M(cont ' d) 

Output (cont’d) 

CMY  Pitching  moment  about  y  axis 

CMZ  Pitching  moment  about  z  axis 

DXN  Center  of  pressure  location 

CM25  Pitching  moment  about  quarter  chord  of  MAC 

CMOO  Pitching  moment  about  leading  edge  of  MAC 

Subroutine  Called:  None 

Error  Returns:  None 
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Subroutine  INFLU 


Purpose: 

Method : 


Use: 


To  calculate  the  three  components  of  velocity  induced  by 
a  constant  source  distribution  on  a  given  panel. 

The  method  of  Hess  and  Smith  (Reference  1)  is  used.  The 
velocity  component  formulas  are  summarized  in  the  Aero¬ 
dynamic  Theory  Section  of  this  report. 

Call  INFLU  (I,  J,  XE,  YE,  XN,  YN,  ZN,  XC,YC,  AC,  A,  T, 

F,  VX,VY,  VZ,  SYM). 

Input : 

I  Control  point  number 

J  Influencing  panel  number 

XE  x  coordinate  of  panel  corner  point  in  panel 

coordinate  system 

YE  y  coordinate  of  panel  corner  point  in  panel 

coordinate  system 

XN  x  coordinate  of  control  point  I  in  reference 

coordinate  system 

YN  y  coordinate  of  control  point  I  in  reference 

coordinate  system 

ZN  z  coordinate  of  control  point  I  in  reference 

coordinate  system 

XC  x  coordinate  of  centroid  of  panel  J  in  reference 

coordinate  system 

YC  y  coordinate  of  centroid  of  panel  J  in  reference 

coordinate  system 

ZC  z  coordinate  of  centroid  of  panel  J  in  reference 

coordinate  system 

A  Transformation  matrix  of  panel  J 

T  Maximum  diagonal  of  panel  J 

F  Area  of  panel  J 
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Subroutine  INFLU (cont ’d) 

Input : (cont  *d) 

SYM  Logical  variable  denoting  symmetry  in  panels  about 

x,  z  plane.  (SYM  is  true  if  the  panel  and  control 
point  lie  on  the  same  side  of  the  plane  of  symmetry; 
SYM  is  false  otherwise.) 


Error  Returns:  Calls  EXIT  if  control  point  lies  on  edge  of 

panel. 
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Subroutine  LATICE 


Purpose: 

Method: 

Use: 


To  calculate  the  three  components  of  velocity  induced  by  a 
given  vortex  lattice. 

The  method  of  Rubbert  and  Saaris  (Reference  2)  is  used. 

The  velocity  component  formulas  are  summarized  in  the 
Aerodynamic  Theory  section  of  this  report. 

Call  LATICE  (XGI,  YGI,  ZGI,  XC2,  ZG2 ,  GA,  N,  XP,  YP,  ZP, 

U,  V,  W,  SYM,  IL) 

Note:  All  points  are  given  in  the  reference  coordinate  system 
Input : 

N  Number  of  bound  vortices  in  lattice 

XGI  x  coordinate  of  inboard  end  of  bound  vortex 

YGI  y  coordinate  of  inboard  end  of  bound  vortex 

ZGI  z  coordinate  of  inboard  end  of  bound  vortex 

XG2  x  coordinate  of  outboard  end  of  bound  vortex 

YG2  y  coordinate  of  outboard  end  of  bound  vo>-' 

ZG2  z  coordinate  of  outboard  end  of  bound 

GA  Relative  strengths  of  bound  vortices  in  iaLi-  'C  ' 

XP  x  coordinate  of  control  point 

YP  y  coordinate  of  control  point 

ZP  z  coordinate  of  control  point 

SYM  Symmetry  parameter 

IL  Vortex  Number 

Output : 

U  x  component  of  induced  velocity  at  control  point 

(in  reference  coordinate  system) 

V  y  component  of  induced  velocity  at  control  point 

(in  reference  coordinate  system) 
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Subroutine  LATICE(cont ' d) 

Output: (cont 'd) 

W  z  component  of  induced  velocity  at  control  point 

(in  reference  coordinate  system) 

Subroutine  Called:  None 

Error  Returns:  Calls  EXIT  if  control  point  lies  on  a  vortex 

line. 
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Subroutine  LINCOF 


Purpose: 

Method: 

Use: 


Subroutines 


To  compute  the  equation  for  the  given  line  segment. 

The  subroutine  LINCOF  computes  the  coef ficientes  in  the 
equation  for  the  given  line  segment.  The  equation  used 
is: 


AV  +  BW  +  C  «  0 


Call  LINCOF  (V^V^.V^W^  A,  B, 

C) 

Input : 

vr  wi 

The  coordinates  of  the 

first  point 

v2,  w2 

The  coordinates  of  the 

second  point 

Output: 

A.B.C, 

The  coefficients  in  the  equation 

Called: 

None 

Error  Returns: 


None 


Subroutine  OVRLAP 


Purpose:  To  eliminate  overlapped  panels  and  overlapped  parts  of 

panels. 

Method:  The  subroutine  OVRLAP  goes  through  two  phases  in  order  to 

eliminate  the  hidden  lines.  The  first  phase  runs  through 
each  panel  listed  in  the  array  LL  and  tests  to  see  if  it 
overlaps  or  is  overlapped  by  any  of  the  subsequent  panels 
in  the  list  LL.  If  one  of  the  panels  is  completely  over¬ 
lapped,  it  is  eliminated.  If  one  of  the  panels  is  partially 
overlapped  by  the  other,  the  overlapping  panel  is  placed  on 
a  list  (L)  of  such  panels  associated  with  the  overlapped 
panel . 

The  second  phase  runs  through  the  list,  L,  of  partially 
overlapping  panels,  and  plots  the  nonhidden  line  segments. 
Once  their  lines  have  been  plotted  the  panel  is  eliminated 
from  the  list  LL. 

Use:  Call  OVRLAP (NN,  LNEXT) 

Input : 

NN  The  number  of  vertices 

LNEXT  The  number  of  panels  in  LL 

In  addition  to  the  parameters  in  the  calling  sequence, 
the  subroutine  OVRLAP  uses  the  data  in  the  common  arrays 
V,  W,  DIST,  and  LL. 

Output: 

OVRLAP  draws  all  the  partially  overlapped  panels,  and  also 
modifies  the  array  LL. 

Subroutines  Called:  LINCOF  Used  to  compute  the  equation  for  a  line 

segment 

PLOT  Used  to  draw  the  nonhidden  line  segments 

Error  Returns:  There  are  four  error  stops  in  the  subroutine 

OVRLAP 

Stop  4 
Stop  5 
Stop  10 
Stop  12 


91 


Subroutine  PLTPAN 


Purpose:  To  draw  the  panels. 

Method:  The  subroutine  PLTPAN  has  three  modes  of  operation.  The 

first  mode  (10P  ■  0)  simply  initializes  the  array  IVECT. 
The  second  mode  (10P  *  2  or  3)  enters  a  sequence  of 
connected  line  segments  into  the  array  IVECT  (10P  *»  3 
for  the  starting  point  of  the  first  segment,  and  10P  *  2 
for  the  subsequent  points) .  The  third  mode  (10P  =  1) 
loops  through  the  array  IVECT  and  draws  the  indicated 
line  segments. 

Use:  Call  PLTPAN  (L,  10P) 

Input: 

L  The  index  of  the  coordinates  of  the  point  in 

question 

10P  The  mode  parameter  (see  above) 

The  data  in  the  common  arrays  V  and  W  are  used  in  addition 
to  the  variables  in  the  calling  sequence 

Output : 

The  output  of  this  routine  is  the  plot 
Subroutines  Called:  PLOT  Used  to  draw  the  panels 

Error  Returns:  None 
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Subroutine  SLEQW 


Purpose: 


To  solve  a  system  of  linear  equations  by  direct  inversion. 


Method:  Gaussian  algorithm  for  solution  of  a  system  of  linear 

equations  with  pivoting. 

Use:  Call  SLEQW  (A,  MM,  R,  MN,  M,  N,  ILL) 

Input : 

A  Matrix  of  coefficients  of  equations  (dimensioned 

MM  x  MM  in  calling  program) 

R  Right  side  vector  matrix  (dimensioned  MM  x  MN  in 

calling  program) 

MM  Maximum  dimensions  of  A 

MN  Maximum  number  of  right  side  vectors 

M  Actual  dimensions  of  A 

N  Number  of  right  side  vectors 


Output : 


Solution  vector 


Subroutine  Called: 


Error  Returns: 


ILL  *  -1  if  system  of  equations  is  ill  conditioned 
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Subroutine  SOLVE 


Purpose: 

Method : 

Use: 


To  solve  a  system  of  linear  equations  by  an  iterative 
procedure. 

The  system  of  linear  equations  is  solved  using  the  Gauss- 
Seidel  iterative  procedure,  with  direct  solution  of  the 
vortex  lattice  partition.  The  method  is  described  in 
the  section  of  this  report  titled,  "The  Boundary  Condition 
Equations"  (see  page  9). 

Call  SOLVE  (A,  B,  X,  HA,  HB,  NS,  N,  LHA,  F,  EPS,  1W,  NIT, 
TP 10,  ITA,  ILL,  HH) 

Input : 

A  Row  of  influence  coefficient  matrix 

B  Right  side  of  boundary  condition  equation 

X  Array  of  source  and  vortex  strengths  (solution 

input  vector) 

HA  Vortex  lattice  influence  coefficient  matrix 

HB  Right  side  of  vortex  lattice  equations 

NS  Number  of  source  panels 

N  Number  of  source  panels  and  vortex  lattices 

(Maximum  size  of  matrix  A) 

LHA  Maximum  dimension  of  matrix  HA 

F  Relaxation  factor  (set  equal  to  unity) 

EPS  Solution  residual  limit 

IW  Initial  value  switch 

If  IW  =  1,  X(I)  =  0 

If  IW  =  0,  X(I)  obtained  from  previous  solution 

NIT  Maximum  number  of  iterations 

TPIO  Name  of  file  used  for  storing  matrix  A 

ITA  Number  of  x  values  printed  out  if  ILL  =  2.  If 

ITA  =  0,  all  x  values  are  printed. 
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Subroutine  Called:  SLEQW 

Error  Returns:  If  ILL  *  1,  subroutine  writes  error  message,  and 

returns.  The  subroutine  calls  EXIT  if  the  solution 
diverges. 
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Program  WBAERO 


Purpose: 


Method : 


Use: 


To  calculate  the  pressure  coefficients  at  the  panel  control 
points  on  wings,  bodies,  and  wing-body  combinations  in 
subsonic  compressible  flow. 

The  panel  corner  points  computed  by  subroutine  WBPAN  are 
read  from  the  auxiliary  file  TAPE  11.  Subroutine  ELEMEN 
is  then  called  for  each  panel  in  turn.  It  calculates  the 
control  point,  the  transformation  matrix,  and  transforms 
the  corner  points  from  the  reference  coordinate  system 
to  the  panel  coordinate  system.  The  panel  control  points 
and  corner  points  are  then  transformed  to  the  analog  body 
using  subroutine  ANALOG,  in  preparation  for  calculation 
of  the  aerodynamic  influence  coefficients.  Subroutine 
INFLU  calculates  the  three  components  of  velocity  induced 
by  the  source  panels,  and  subroutine  LATICE  calculates 
the  three  components  of  velocity  induced  by  each  vortex 
lattice.  These  velocities  are  combined  to  form  the  matrix 
of  aerodynamic  influence  coefficients,  one  row  at  a  time. 
The  influence  coefficient  matrix  is  stored  on  auxiliary 
file  TAPE  10  in  row  order,  and  the  three  components  of 
velocity  are  stored  on  TAPE  12,  also  in  row  order.  The 
right  side  of  the  boundary  condition  equation  is  computed 
for  each  angle  of  attack  and  yaw,  and  the  system  of 
equations  solved  for  the  source  and  vortex  strengths  by 
calling  subroutine  SOLVE.  A  detailed  description  of  the 
method  is  given  in  the  Aerodynamic  Theory  section  of  this 
report.  The  pressure  coefficients  are  then  obtained  by 
summing  the  products  of  the  velocity  components  and 
singularity  strengths,  and  applying  equations  (26)  and  (27) 
given  in  the  Aerodynamic  Theory  section. 

Call  OVERLAY  (FRWB,  3,  0) 

Input : 

TEXT  Identifying  title 

PRINT  Print  option  selected  (see  input  section) 

NIT  Maximum  number  of  iterations 

IEPS  Exponent  of  10  setting  limit  on  residue  of 
iterative  solution 

ITYPE  Type  of  solution  procedure  selected 
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Program  WBAERO(cont 'd) 


Input : (cont ’d) 

ISAVE  Control  parameter  for  reading  influence  coefficients 
and  velocities  components  for  auxiliary  files 
TAPE  10  and  TAPE  12,  on  to  TAPE  11. 

SIM  Panel  symmetry  parameter 

KUT  Vortex  lattice  control  parameter 

NBV  Number  of  body  vortices  having  same  strength  as 

adjacent  wing  vortices 

NV  Number  of  wing  vortices  associated  with  each  body 

vortex 

KOMPR  Compressibility  rule  parameter 

POINTS  Number  of  field  points  requested 

NORPAN  Number  of  panels  with  non-zero  normal  velocity 

NMA  Number  of  Mach  numbers 

MA  Mach  number 

NAL  Number  of  angles  of  attack  or  yaw 

ALPHA  Angle  of  attack 

BETA  Angle  of  yaw 

X,Y,Z  Source  panel  corner  points  in  reference  coordinate 
system 

XL,YL,  Vortex  panel  corner  points  in  reference  coordinate 
ZL  system 

XP,YP,  Field  point  coordinates 

ZP 

GA  Relative  strengths  of  bound  vortices  in  vortex 

lattices 


J 


Panel  Number 
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Program  WBAERO(cont ' 
Output: 
LP 

BETAl 

NP 

NORVEL 

XC, YC, 
ZC.XCI, 
YCI.ZCI 

XE,  YE 

A 

F 

T 

VX, VY, 
VZ 

AM 

RULE  1 

RULE  2 

VXU.VYU 

VZU 

RS 

ILL 

SIGMA 

VXR.VYR 

VZR 


d) 

(cont'd) 

Bound  vortex  number 
Prandtl-Glauert  factor 

Panel  number  at  which  non-zero  normal  specified 
Normal  velocity 

Panel  control  points  in  reference  coordinate 
system 

Panel  control  points  in  reference  coordinate 
system 

Panel  transformation  matrix 
Panel  area 

Maximum  diagonal  of  panel 

Three  components  of  Induced  velocity  in  reference 
coordinate  system 

Component  of  velocity,  normal  to  plane  of  panel 
(influence  coefficient) 

Gothert's  rule  1  selected 

Gothert's  rule  2  selected 

Three  components  of  the  free-stream  velocity 
vector 

Right  side  of  boundary  condition  equations 
Matrix  solution  Indicator 
Array  of  source  and  vortex  strengths 
Resultant  velocity  component  arrays 
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V  Magnitude  of  resultant  velocity  vector 

CP  Pressure  coefficient  array 

Subroutines  Called:  ELEMEN 

ANALOG 

INFLU 

LATIC 

SOLVE 

FORMOM 

EXIT 

Error  Returns:  Program  calls  EXIT  if: 

1.  MA  >  1.0 

2.  J  >  1500 

3.  LP  >  40 

4.  IL  >  35 


5.  ILL  -  1 
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Program  WBOLAY 

Purpose:  Main  overlay  for  wing-body  analysis  program 

Method:  To  call  the  primary  overlay  programs  WBPAN  and  WBAERO 

Use:  OVERLAY  (FRWB,  0,  0) 

FRWB  is  overlay  file  name 

Program  Called: 

OVERLAY  (FRWB,  1,  0)  (WBPAN) 

OVERLAY  (FRWB,  2,  0)  (WBPLOT) 

OVERLAY  (FRWB,  3,  0)  (WBAERO) 

Subroutines  Called:  Exit 
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Program  WBPAN 


Purpose: 

Method : 


To  calculate  panel  subdivision  for  wings,  bodies,  or 
wing-body  combinations. 

Several  alternate  paths  are  available  in  this  subroutine 
depending  on  the  values  of  the  control  parameters  selected. 
Individual  panel  corner  point  arrays  are  read  in  if  SINGPA  ■ 
1,  body  section  data  is  read  in  if  CASE  ■  1  or  3,  and  wing 

section  data  is  read  in  if  CASE  ■  2  or  3. 

Four  options  are  available  for  reading  in  the  body  section 
area.  If  OPT  ■  0,  the  y  and  z  coordinates  of  the  panel 
corner  points  are  required.  If  OPT  -  1,  the  preceding 
section  comer  points  are  used,  and  no  additional  data 
is  read.  If  OPT  ■  2,  the  panel  corner  points  are  specified 
by  the  polar  coordinates  r  and  0.  For  bodies  of  revolution 
having  uniform  panel  spacing,  OPT  *  3  provides  a  simplified 
input  option,  and  requires  the  input  of  only  the  radius  and 
8  increment  for  each  section.  If  OPT  ■  0  or  2  have  been 

selected,  the  x  coordinate  of  the  panel  corner  points  may 

be  shifted  out  of  the  plane  of  the  section  by  an  amount 
Ax  to  allow  more  freedom  in  paneling  intersections.  In  all 
cases,  the  program  calculates  the  x,y,z  coordinates  of  the 
four  panel  comers  in  the  reference  coordinate  system,  and 
writes  them  on  the  auxiliary  file  TAPE  11  and  the  output 
file. 

The  wing  section  data  is  input  as  airfoil  coordinate  arrays. 
These  arrays  may  be  given  in  the  reference  coordinate  system, 
or  in  terms  of  the  local  percent  chord.  If  the  latter  option 
is  selected,  the  chord  length,  twist  angle,  and  twist  center 
must  be  specified  for  each  section.  An  arbitrary  dihedral 
angle  may  also  be  specified  for  each  section. 

The  internal  vortex  lattice  panels  are  located  on  the  mean 
camber  line  of  the  wing  section.  The  relative  strength 
(GAMMA)  array  must  be  specified  for  each  section.  In 
addition,  the  y  coordinate  of  panel  corner  points  may  be 
shifted  out  of  the  plane  of  the  section  by  an  amount  Ay 
to  allow  more  freedom  in  paneling  wing  tips  and  wing- 
body  intersections.  In  all  cases,  the  program  calculates 
the  x,y,z  coordinates  of  the  four  corners  of  the  surface 
panels  and  vortex  lattice  panels  in  the  reference  coordinate 
system,  and  writes  them  on  the  auxiliary  file  TAPE  11 
and  the  output  file.  For  wing-body  combinations  (CASE  ■  3) 
additional  vortex  lattice  panels  may  be  specified  inside  the 
body.  The  input  required  is  similar  to  that  described  above 
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Program  WBPAN(cont '  d) 
Method: (cont ' d) 


Use: 


the  wing  panels.  The  program  calculates  the  x,y,z 
coordinates  of  the  four  corners  of  the  additional  vortex 
panels  in  the  reference  coordinate  system,  and  writes  them 
on  the  auxiliary  file  TAPE  11  and  the  output  file. 

Call  OVERLAY  (FRWB,  1,  0) 


TEXT 

CASE 

PLOT 

SIM 

I  SAVE 

PRINT 

SINGPA 

NOPAN 


XX, YY, 
ZZ 

NB 

NW 

NV 

XBE.YBE 

ZBE 

MB,NW 


MV 

OPT 

CHD 


Identifying  title 
Component  identification  parametet 
Plot  selection  parameter 
Configuration  symmetry  parameter 
Save  tape  parameter 
Print  option 

Single  panel  selection  parameter 

Number  of  panels  to  be  deleted 

Panel  corner  point  coordinate  in  reference 
coordinate  system 

Number  of  body  sections 

Number  of  wing  sections 

Number  of  vortex  lattices  in  body 

Section  coordinates  in  reference  coordinate 
system 

Number  of  panel  corner  points  in  wing  or  body 
section 

Number  of  vortex  panels  in  body  vortex  lattices 
Corner  point  input  option  (see  description  above) 
Panel  chord 
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Program  WBPAN(cont ' d) 


In:  ut : 

ALF 

XAL 

KOORD 

FLAG 

PEL 

DELTA 

YO.ZO 

THET 

A 


B 


C 

D 

WAKE 


(cont 'd  ) 

Panel  twist  angle (degrees) 

Twist  reference  point 

Wing  panel  coordinate  system  parameter 

Branch  point  indicator (required  to  change 
number  of  panels  from  section  to  section) 

Dihedral  parameter 

Wing  section  dihedral  (degrees) 

Coordinate  of  axis  of  rotation  for  dihedral 

Theta  increment  (degrees)  for  OPT  ■  3 

-z  coordinate  of  body  panel  corner  point  if 
OPT  =  0  and  CASE  >=  1 

-r  coordinate  of  body  panel  corner  point  if  OPT  >  0 
and  CASE  *  1 

-z  coordinate  of  wing  panel  corner  point  if 
CASE  -  2 

-y  coordinate  of  body  panel  corner  point  if 
OPT  =  0  and  CASE  =  1 

-8  coordinate  (in  degrees)  of  body  panel  corner 
point  if  OPT  >  0  and  CASE  ■  1 

-x  coordinate  of  wing  panel  corner  point  if  CASE  ■ 

2 

-Relative  strength  of  vortex  lattice  panels  (GAMMA) 

-Ax  shift  of  body  panel  corner  point  if  CASE  -  1 

-Ay  shift  of  wing  panel  corner  point  if  CASE  *  2 

Length  of  vortex  lattice  in  wake  in  percent  of 
local  chord 
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Program  WBPAN(cont 'd) 


Input : 
POINT 

XLP,YH 

ZLP 

Output 

PANEL 

NPAN 

XX,  YY 

ZZ,X,Y 

Z 

XH,YH, 

ZH.XC, 

YC.ZC 

MINUS  : 

Subroutine  Called: 
Error  Returns: 


(cont  *d) 

Location  of  vortex  lattice  control  point  behind 
trailing  edge  in  percent  of  local  chord 

,  Coordinates  of  terminal  points  of  streamwise 
vortices  (input  only  if  FLAG  =  3  and  CASE  >  2) 


Panel  number 

Panel  numbers  of  panels  to  be  deleted 

Panel  corner  point  coordinates  in  reference 
coordinate  system 

Vortex  lattice  panel  coordinates  in  reference 
coordinate  system 

End  of  record  mark  for  TAPE  11 
EXIT 

Program  calls  EXIT  if: 

1.  NB  >  70 

2.  MB  >  60 

3.  NW  >  AO 
A.  NW  >  60 

5.  NV  >  AO 

6.  MV  >  60 
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Program  WBPLOT 


Purpose:  To  plot  the  panel  geometry 


Method:  The  plot  parameters  and  viewpoint  coordinates  are  read 

from  the  input  file,  and  the  panel  corner  point  coordi¬ 
nates  are  read  from  TAPE  11.  The  data  is  stored  in 
labelled  COMMONS  SCRAT,  PLODAT,  and  PLOPAR  prior  to 
calling  the  plot  subroutines. 


Use:  Call  OVERLAY  (FRWB,  2,0) 

Input : 

NVU  Number  of  viewpoints 

IPRINT  Print  parameter 

IHIDE  Hidden  line  parameter 

TBUG  Debug  print  parameter 

VUE  Viewpoint  coordinate 

X,Y,Z  Panel  corner  point  coordinates 


Output: 

M  Number  of  panels 

N  Number  of  corner  points 

L  Array  of  corner  point  indices 


Subroutines  Called:  PLOTS 

FACTOR 

DRAW3D 

PLOT 


Error  Returns: 


1.  NVU  >  4  -  Subroutine  writes  error  message  & 

returns 

2.  M  >  1500  -  Subroutine  writes  error  message  & 

returns 

3.  N  >  3000  -  Subroutine  writes  error  message  & 

returns 
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APPENDIX  III 


SAMPLE  INPUT 


The  input  for  a  wing-body-tail  configuration  is  given  in  this  appendix. 

Figure  28  is  a  computer  generated  plot  of  the  paneling  used  for  this 
configuration.  Panels  are  input  for  one  side  of  the  configuration 
only . 

The  body  is  an  ogive  cylinder  with  a  blunt  base,  and  is  subdivided  into 
72  panels,  including  base  panels.  The  wing  has  a  truncated  delta  plat¬ 
form  and  a  thickness/chord  ratio  of  approximately  15  percent.  It  is 
subdivided  into  69  panels.  The  vertical  tail  has  a  swept  constant  chord 
planform  and  the  same  airfoil  section  as  the  wing.  It  contains  an 
additional  32  panels.  The  total  number  of  source  panels  on  the  configu¬ 
ration  is  168. 

Vortex  lattice  panels  are  automatically  included  inside  the  wing,  while 
the  vortex  lattices  inside  the  vertical  tail  are  input  as  body  vortices. 
Eight  vortex  lattices  are  used  in  the  wing,  body,  and  vertical  tail. 

The  sample  input  is  chosen  to  illustrate  some  of  the  special  features 
of  the  program.  These  features  are  described  below: 

1.  Plot  Option  -  Tha  plot  option  is  selected  by  setting  PLOT  ■  1 
on  Card  2.  The  plot  parameters  are  set  on  Cards  9  and  10. 

2.  Body  Input  -  A  simplified  body  input  definition  is  obtained 
by  setting  OPT  *  3  on  Card  4B.  This  option  only  applies  to 
bodies  of  revolution. 

3.  Wing  and  Tail  Input  -  The  wing  input  card  set  includes  the 
vertical  tail.  The  main  wing  uses  the  standard  input  option, 
while  the  vertical  tall  is  defined  in  a  horizontal  plane  and 
rotated  into  the  vertical  position  using  the  dihedral  option 
(DEL  ■  1  on  Card  5B) .  Since  the  wing  and  tail  have  the  same 
airfoil  section,  this  data  is  input  only  for  the  first  wing 
section,  and  the  remainder  are  defined  by  setting  OPT  ■  1  on 
Card  5B.  The  vortex  lattices  are  automatically  calculated 
for  the  main  wing,  but  omitted  from  the  vertical  tail  by 
setting  FLAG  =  2.  The  vertical  tail  vortex  lattices  are 
input  later  as  body  vortices.  The  trailing  vortices  from 
the  wing  vortex  lattice  are  extended  100  units  into  the 

wake  using  Card  Set  7,  and  setting  FLAG  ■  3  on  Card  5B. 

The  vortex  lattice  control  points  are  defined  by  Card  6. 
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4.  Body  Vortices  -  A  vortex  lattice  is  input  inside  the  body  to 
provide  a  mechanism  for  carry-over  of  lift.  The  strengths  and 
locations  of  the  bound  vortices  in  this  lattice  are  chosen  to 
match  those  of  the  adjacent  wing  vortices  in  this  lattice,  and 
the  trailing  vortex  is  extended  into  the  wake  such  that  it 
terminates  at  the  same  location  as  the  inboard  wing  trailing 
vortices. 

Body  vortices  are  also  input  inside  the  vertical  tail.  The  bound 
vortices  are  located  in  the  plane  of  symmetry  under  the  spanwise 
panel  edges,  and  given  a  strength  proportional  to  the  airfoil  local 
thickness  distribution.  The  trailing  vortices  are  extended  100 
units  into  the  wake,  using  Card  Set  8D-1.  The  vertical  tail  vortex 
lattice  control  panels  are  defined  by  setting  OPT  B  2  on  Card  8B, 
and  reading  in  the  control  point  coordinates  on  Card  8D-2  and  3. 
Finally,  the  symmetry  option  is  ignored  for  these  vortex  lattices 
by  setting  SIMOPT  *  1  on  Card  8B.  This  suppresses  the  image  vortex 
lattice  system  which  in  this  case  would  exactly  cancel  the  defined 
vortex  lattice  system  in  the  vertical  tail  and  result  in  a  singular 
matrix  being  generated  in  the  aerodynamic  section  of  the  program. 

5.  Field  Point  uption  -  The  field  point  option  is  selected  by  setting 
POINTS  =  16  on  Card  16.  The  coordinates  of  the  16  field  points 
follow  on  Card  Set  16A. 

6.  Normal  Velocity  Option  -  The  normal  velocity  option  is  selected  by 
setting  NORPAN  =  8  on  Card  16.  The  indices  of  the  8  points  are 
identified  on  Card  Set  16B.  In  this  example,  NORVEL  -  0,  so  the 
normal  component  of  the  onset  velocity  is  set  equal  to  zero  on 

the  blunt  base  of  the  body. 
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APPENDIX  IV 


SAMPLE  OUTPUT 


The  output  for  a  wing-body-tail  configuration  is  given  in  this  appendix. 
Only  standard  output  is  given,  since  the  print  option  PRINT  *  0  was 
selected. 


1.  List  of  Input  Data  -  The  program  lists  all  input  data. 

2.  Panel  Corner  Point  Coordinates  -  The  corner  point  coordinates 
of  the  body  panels,  wing  panels,  and  vortex  panels  are  listed. 

3.  Aerodynamic  Parameters  -  Selected  aerodynamic  parameters  are 
listed,  together  with  the  CPU  time  required  to  calculate  the 
aerodynamic  matrix  and  solve  the  system  of  equations. 

4.  Solution  Parameters  and  Singularity  Strengths  -  A  list  giving 
the  residuals  obtained  after  each  step  of  the  iteration 
procedure  is  printed,  together  with  the  final  array  of 
singularity  strengths  obtained,  in  order  of  panel  number. 

5.  Velocity  and  Pressure  Distribution  -  The  three  components  of 
velocity,  the  magnitude  of  the  velocity  vector,  and  the 
pressure  coefficient  are  tabulated  for  each  panel,  together 
with  the  location  of  the  panel  control  point. 

6.  Total  Coefficients  -  The  axial,  normal,  and  side  forces,  the 
moments  about  the  three  coordinate  axes,  the  pitching  moment 
about  the  reference  point  and  the  quarter  chord  of  the  MAC, 
the  center  of  pressure,  and  the  lift  and  drag  coefficients 
are  listed. 

Items  4,  5,  and  6  are  repeated  for  each  angle  of  attack  or  yaw 

selected. 
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LIST  OF  SYMBOLS 


A  Aerodynamic  matrix 

a  Aerodynamic  influence  coefficient 

B  Prandtl-Glauert  factor 

C  Aerodynamic  coefficient 

c  Chord 

D  Perpendicular  distance  from  panel  edge  to  control  point 

d  Distance  between  panel  corner  points 

J  Geometrical  parameter  for  source  panel 

L  Length  of  line  vortex 

M  Mach  number,  pitching  moment 

N  Number  of  singularities 

NL  Number  of  vortex  lattices 

NS  Number  of  source  lattices 

n  Direction  cosine  of  normal  vector 

P,Q  Geometrical  parameters  for  source  panels 

q  Magnitude  of  resultant  velocity  vector 

R  Component  of  free-stream  velocity  vector  normal  to  panel 

r  Distance  from  panel  corner  point  to  control  point 

S  Area 

T  Geometrical  parameter  for  source  panel 

t  Component  of  transformation  matrix 

U,V,W  Components  of  resultant  velocity  vector 
u,v,w  Perturbation  velocity  components 
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List  of  Symbols  (cont'd) 


V  Velocity  vector  magnitude 

x(y,z  Cartesian  coordinates  of  point 
a  Angle  of  attack 

3  Angle  of  yaw 

y  Vortex  strength 

e  Small  reference  value 

k  Ratio  of  specific  heats  for  air 

0  Angle  between  panel  edge  and  line  joining  panel  corner  to 
control  point 

o  Singularity  strength 

s  Source  strength 

£,n  Panel  coordinates 

Mathematical  Symbols 
-*■  Vector 

j  Absolute  value 

x  Cross  product 

•  Scalar  product 

Average  value 
[  J  Matrix  quantity 

{  }  Vector  array 

I  Summation 

1  Integral 

A  Incremental  value 
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List  of  Symbols  (cont'd) 


Subscripts 

a  Analogous  body 

D  Drag 

i  Panel  control  point 

j  Panel  corner  point 

L  Lift 

M  Moment 

P  Pressure 

s  Source 

S  Side 

v  Vortex 

W  Wing 

X  Axial 

Y  Lateral 

Z  Normal 

x,y,z  Reference  axis  direction 
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